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Abstract. We consider the approximation scheme of the American call op- 
tion via the discrete Morse semiflow. It is the minimizing scheme of a time- 
semidiscretized variational functional. In this paper we obtain a rate of conver- 
gence of approximate solutions. In addition, the convergence of approximate 
free boundaries is proved. 

1 Introduction 

In this paper we consider an approximation scheme to the following obstacle problem: 

mm ^-C T - ^-S 2 C SS ~ (r - q)SC s + rC,C - ^ =0 in (0,T) x (0,+oo), 

C(T, S) = <S>(S) := max(5 - K, 0) for S G [0, +oo), 

C(r,0) = forrG(0,T), 

The above equation is called the Black-Scholes equation for the American call option. Here 
C — C(t, S) is the option price, the positive constants a, r, q, K denote, respectively, the 
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volatility, the interest rate, the dividend and the strike price. Throughout this paper we 
set a = for simplicity and assume r > q. 

It is known that (11. ip has a unique free boundary. Indeed, by van Moerbeke [34] . 
we see that under the assumption r > q, there exists a unique S* G C 1 (0,T) such that 
(S*)' < and 

{S | C(r, 5) > = (-00, S*(t)), {S I C(r, 5) = = [5"(r), +00) 

for each r G (0, T), 

(1.2) limSVr) = — (>K). 

t/t q 

See Wilmott - Howison - Dewynne [551 P-124] for the formal derivation of (II. 2p . The 
family {5 , *(r)} < r <T is called the free boundary or the optimal exercise boundary; for 
each r G (0, T), S*(t) indicates the value of the current stock price under which the 
holder of the option should (optimally) exercise it. 

From the viewpoint of mathematical finance, it would be very convenient to obtain 
{S*(t)} 0<t< t explicitely. However, this seems to be very difficult and thus many people 
have studied numerical schemes for (11.11) . especially to approximate the free boundary. 
Brennan - Schwartz [7] introduced a fully implicit difference scheme for the American 
put option and obtained a numerical solution. The convergence of their scheme was 
proved by Jaillet - Lamberton - Lapeyre [13j in the framework of variational inequali- 
ties. Lamberton [TB] considered the binomial tree method and the finite difference one to 
approximate {S*(t)} 0<t< t and showed the convergence of the approximate free bound- 
ary by the probabilistic argument and the analytical one. He also obtained in [T9l 120] 
some error estimates for the stochastic approximation to the optimal stopping problems 
including the American options. Amin - Khanna |2j treated a discrete time model for 
the American option and proved the convergence of the discrete American option value 
to the continuous one. Jiang - Dai [T4] obtained similar results to those in [T5] by the 
method of viscosity solutions. Omata - Iwasaki - Nakane - Xiong - Sakuma [29] proposed 
an approximation scheme to (11.11) different from the above ones and obtained a numerical 
result. 

The approximation scheme by [2H] is based on the discrete Morse semiflow (DMS), con- 
sisting of the minimization of a time-semidiscretized variational functional. The DMS was 
first used by Rektorys [31] to obtain the solutions of linear parabolic equations. Kikuchi 
[T5| fT6] applied the DMS to construct the solutions of parabolic equations associated with 
a variational functional of a harmonic map type. Besides, in [321 I23~l 1241 125] Nagasawa 
and Tachikawa used the DMS to show the existence and asymptotic behavior of solutions 
of some semilinear hyperbolic systems. Nagasawa - Omata [22J considered the behavior 
of the DMS for a free boundary problem. Some applications of the DMS to numerical 
analysis have been treated in Omata [261 EZ], Omata - Okamura - Nakane [30] and Omata 
- Iwasaki - Kawagoe [28] . 

The purpose of this paper is to discuss the convergence of the approximation scheme 
by [2J5]. Our results are a rate of convergence of the approximate solutions and the 
convergence of the approximate free boundary. The former result is obtained by applying 
the rate of convergence of product formula for semigroups by Bentkus - Paulauskas [6] 
and the precise comparison argument for viscosity solutions by Ishii - Koike [TT] , in which 
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they obtained the rate of convergence in elliptic singular perturbations. The latter one is 
proved by the limit operation of viscosity solutions due to Barles - Perthame [3J H]. 

This paper is organized as follows. In Section 2 we introduce the approximation scheme 
by [29] and state the main results. In Section 3 we discuss the solutions of (12. 2p below. In 
Subsection 3.1 we briefly prove the existence and uniqueness of solutions. In Subsection 
3.2 we show some properties of solutions. Section 4 is devoted to the DMS associated 
with (D (DMS-BS for short) and the free boundary of the DMS-BS. In Subsection 4.1 
we derive some estimates for the DMS-BS. Subsection 4.2 is devoted to the existence and 
uniqueness of the free boundary of the DMS-BS. In Subsection 4.3, we give a proof of 
Theorem 14.51 in subsection 4.1, an estimate of the difference of the DMS-BS. In Section 
5 we prove our main results. Section 6 is the Appendix. In Subsection 6.1 we discuss 
the formal asymptotic expansion of an ODE related to (14. 3p below. This expansion is 
used to construct sub- and super solutions of (12.41) and (I4.3p . In Subsection 6.2 we give an 
estimate for some coefficients appearing in the estimate of Theorem 14.61 

In the following of this paper, we denote by C various constants depending only on 
known ones. The value of C may vary from line to line. 



2 Approximation scheme and Main Results 

In this section we state the approximation scheme by [29] and our main results. 

We reformulate (11.11) in the following way. Put t := T — r, x := \og(S/K), a : = 
(r - q - l)/2 and U(t,x) := S a C(r, S)/K a+1 . Then flU} turns to 

min {U t - U xx +[3U,U-(p} = in (0, T) x R, 
(2.1) { U(0, x) = <p(x) := e ax max(e I - 1, 0) for x ER, 

U(t,x) — >0 (ac->-oo) forte (0,T), 

where j3 := a 2 + r. From the viewpoint of the numerical analysis, we had better restrict 
the problem (12. ip on a bounded interval with respect to x. This restriction seems to be 
reasonable. Because the free boundary {log(S*(T — t)/K)} 0<t<T for (12. ip is bounded for 
each T > and it is easily seen that 

U(t, x) = 0(e JX ) as x — >• — oo for all t G (0, T) and some 7 > 0. 

Hence, putting Q :— (— 1, 1), we consider the following problem instead of (12.11) : 

{min {ut — u xx + /3u, u — ip} = in (0, T) x Q, 
u(0,x)=<p(x) forxeH, 
u{t,±l) = ip{±l) forte(0,T). 

We assume q < r < qe and denote by {x*(t)}o <t< T the free boundary for (12.21) . Note 
by [H] that 



(2.3) 



x* G CHO, T), (x*)'(t) > 0, \imx*(t) = x := log ( - J (G (0, 1)). 

t\o \qj 
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The approximation scheme by |29j is stated as follows. Fix a time step h > 0. Put 
«o := v 9 an d let [r] be the Gauss symbol for rGl For m — 1, 2, ... , [T/h], we consider 
the minimization problem of the following functional: 

J m {u) := ^ J | — ^_±L _|_ ( Ux y _|_ I for u (z fc^ 

K:={v e H\Q) \ v — <f e H%(Q), v>tp a.e. in Q}. 

We observe by the direct method of calculus of variation that there is a unique minimizer 
u m G /C of J m . Moreover, w m satisfies the elliptic variational inequality: 

(2.4) min j Urn ^ m ~ l - u m>xx + (3u m , u m - y?| = in fi, u m (±l) = <p(±l). 

We call the sequence {« m }^J' the DMS-BS. In addition, there is a unique free boundary 
{^m}m=o ^° the DMS-BS, as will be shown in Subsection 4.2 below. 
Under these settings, we define u h (t,x) and x h (t) by 

(2.5) u h {t,x) := u m (x), x h (t) := x m 

for t G [mh, (m + l)h), x G Q and m — 0,1, ... , [T/h]. 
Then our main results are stated as follows. 

Theorem 2.1 Assume q < r < qe. Then for any 5 > 0, there exist K > and h > 
suc/i £/ia£ /or a// /i G (0, ho), 

sup |w fc (£, —u(t, x) | < fCv^l log h\. 

(t,x)e[o,T~8]xn 

Theorem 2.2 Assume q < r < qe and {x*(t)}o <t< T C fi. T/iere /or any 5 > 0, we have 

lim sup \x h {t) - x*(t)\ = 0. 
fo - >Q te[o,T-<5] 

Remark 2.1 (1) The assumptions r < qe and {x*(t)} <i<r C ft are technical ones. Since 
{x*(t)}o <t< T is bounded for each T > (cf. [SI]), replacing Q with a larger interval such 
that {x*(t)}o <t< T C fi, we can show Theorems 12.11 and 12. 2 [ assuming only r > q. 
(2) The | log/i| term appear in Theorem 12. II by some technical reasons. 



3 Solutions of the problem (12.21) 



This section consists of two subsections. In Subsection 3.1, we consider the existence and 
uniqueness of solutions of (12.21) . In Subsection 3.2, we obtain some regularity of solutions 
of (12.21) . To establish the results in these subsections, we use the penalized problem for 

( u\ - u £ xx + (3u £ + Ce(u £ -<p) = in (0, T)_x to, 
(3.1) < u £ (0,x) = (p{x) for i6fl, 

{ u £ {t, ±1) = cp(±l) for t G (0, T), 

where e > 0, ( E {r) := ({r/e) and ( is a smooth function such that 

C > 0, C" < on M, 
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3.1 Existence and uniqueness of solutions 

In this subsection we prove the following theorem. 

Theorem 3.1 There exists a unique solution u of l\2.^i in the a.e. sense and in the sense 
of viscosity solutions such that u G H /1 ' 2,2 ((0, T) x Q) n C([0,T) x H). 

Remark 3.1 See [H] or [T7| for the definition and the theory of viscosity solutions. 

By [21 Chapter 2, 3], there is a unique solution u £ of (13. ip in the sense that u £ — <p G 
L 2 (0,T; uf G L 2 (0,T;L 2 (ft)) and 

/ {wf0 + u e x (j) x + /3w e + C £ (m £ - <^)0}c?x = for all G #q(0) and a.e. t G (0, T). 
Jn 

In addition, w e G L 2 (0, T; # 2 (ft)) n C([0, T); (7(H)). 

We derive some estimates of « e . By the maximum principle, we get 

(3-2) sup ||w e ||c([o,T)xfI) < IMIc(O)- 

£>0 

By the same arguments as in [5J Chapter 2, Section 2.4] we have 

(3.3) sup(||u e || L oo (0iT;H i (n)) + |K|| L 2 (0iT . L 2 (n)) ) < +oo. 

£>0 

To estimate ( £ (u £ — (p) and u £ xx , we need the following lemma. 

Lemma 3.1 There is M 1 > such that u £ > if - M\e in [0, T) x Cl for all £ > 0. 

Proof. Since y> is Lipschitz continuous and convex in Q, we can show that 

(3.4) / (p x (x)(j) x (x)dx < for all G Hq(Q) satisfying > in f2. 
Jn 

Set u(t,x) := v 9 ( a; ) ~~ Mie. Then, we use the above inequality to obtain 

/ {u t <P + Ux&x + Pvrf> + Ce{u - p)<f)}dx < {(3(f + C(-Mx)}(f)dx 
Jn Jn 

for any t G (0, T) and G -f^o(^) satisfying > in Q. Taking Mi > such that 
/3|Mlc(I}) + C( — Mi) < 0, we easily see that u is a weak subsolution of (13.11) . Hence we 
have u < u £ in [0,T) x Q by the maximum principle. Therefore we obtain the result. □ 

Hence from (I3.3P and Lemma [3.11 we get 

(3-5) sup \\( £ (u £ - <f)\\ c([0tT)x m < +oo, sup \\u £ xx \\ L2{(0;T)xn) < +oo. 

e>0 £>0 

Proof of Theorem 13.11 By ( 13. 2p . (13.31) . (13. 5p and Sobolev imbedding, we can extract a 
subsequence {£ n }t=n £ n \ such that for any X" G (0, T) and A G (0, 1/2), asn-> +oo, 

(3.6) u £n — ► u in C A/2 ' A ([0,T'] x H), 

(3.7) «\ <"J — + a,, u xx ) weakly in (L 2 ((0, T) x ft)) 3 . 

We can see that u is a unique solution of (12.21) in the a.e. sense and in the viscosity sense 
(cf. |5J Chapter 3], [8j and [II])- Thus we complete the proof. □ 
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3.2 Some properties of solutions 

The main results of this subsection are stated as follows. Let u be the solution of (12.21) . 

Theorem 3.2 Assume q < r < qe and let x be given in (\2.3\) . Then u G W^ ,oo ((0, T) x 
Q) and it satisfies the following estimates. 

(1) For any small x\ > 0, there is L\ > such that 

(3.8) |K(*,0IU oo (-a!i^i) < IM*)0IU~(n\(-* 1>a !i)) < In for a.e. t G (0,T), 

(3.9) \\u xx (t, -)|| < llw^^, -)IU°°(fi\(-zi,zi)) < ^i for a.e.t e (0,T). 

(2) There is L 2 > swc/i £/ja£ 

-u(s,y)| < L 2 (|t- s| 1/2 + |x - ?/|) for all (t,x), (s,y) G [0,T) x H. 

Theorem 3.3 T/ie M t nonnegative and continuous in (0, T) x Jl. 

Remark 3.2 (1) In [13], similar results to Theorem 13.21 are obtained in the case Q = Mr. 

(2) Theorem 13.31 is similar to [331 Lemma 5] and 0, Corollary 4.2]. As seen in Section 
5, it plays an important role to prove Theorem 12. 1L 

We prepare some pointwise estimates of solutions of (13.1 j) to prove Theorem 13.21 Let 
u £ be the solution of (13.11) . 

Proposition 3.1 We obtain 

sup |K|U<x. ( ( 0iT)xn) < +oo. 

£>0 

Proof. The barrier construction argument yields that \u e x (t, ±1)| < ||<io||z°°(n) for all 
t G [0, T) and e > 0. We obtain the result by combining the comparison argument for 
viscosity solutions (cf. [121 Section 7]) with Lipschitz continuity of <p and this estimate. 
□ 

Lemma 3.2 Assume q < r < qe and let Xq be defined by <M) . Then there exist x 2 G 
(0, x ) such that u £ (t, x) > <p(x) for all t G (0, T), x G (— 1, x 2 ) and e > 0. 

We can formally show this lemma, according to [351 p. 124]. Let x be defined in (12.31) and 
assume x G (0, 1). It is seen that for small e > 0, 

u\ (t, x) — > Wj(0, x) = f xx (x) — (3tp(x) for all x G (0, 1) as t \ 0. 

It follows from the definition of (p that 

{> if x G (0,x ), 
= itx = x , 
< if x G (x , 1], 
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Thus Lemma [3.21 formally holds with X\ — xq. 

Proof of Lemma 13.21 Step 1. Let u be the classical solution of 

{Ik ~~ ILxx + Pu = in (0, +oo) x Q, 
u(0,x) = (p(x) forxett, 
u(t, ±1) = <p(±l) for t E (0, +oo). 

We prove that 

(3.12) u e (t,x) > u(t,x) > ip(x) for all t > 0, x E (—1, 0] and e > 0. 

Since u is a subsolution of fl3.ll) . it follows from the maximum principle that u < u £ 
in [0, T) x a Hence we get (I3321 by u > in (0, +oo) x and ip = on [-1, 0]. 
Step 2. We show that there exist ti, M 2 > such that 

pt — fig — a; /4s 

(3.13) u(t,x) > <p(x) + M 2 t + / 6 ds for all t E (0, ii] and x E [0, 2xo/3]. 

JO V87TS 



For a > 0, define 

+0O 



(3.14) E a (t, X, y) := ^= ^ { e -(*-y+*™f /4t _ e -(x+y+4an+2a)*/4ty^ 



(3.15) E (t,x,y) := — Le"^ 2 / 44 , z, y) := £7°(t, rr, - E (t,x,y). 

Put E = E 1 and i?i = £^ for simplicity. Then u is given by 

u(t,x) = e~ m [ E(t,x,y)p(y)dy- <p(l) [ e'^ l ~ s) Ey{t - s, x, l)ds. 
Jo Jo 

Differentiating this formula with respect to t, we have 

u t (t, x) = e~ pt [ E t (t, x, y)(f(y)dy - (3e~ pt [ E(t, x, y)^{y)dy - e- pt <f(l)E y (t, x, 1). 
Jo Jo 

We use the facts E t = E xx = E yy , <£>(0) = and the integration by parts to obtain 



u t (t,x) = e~ pt E(t,x,y)& yy (y)-P(p(y))dy + e- f *E(t,x,-W)(p v (-+Q) 
Jo 

= ■■ h,l + h,2- 

We estimate the right-hand side (RHS for short) of the above formula to have (13.131) . 
Some calculations yield that for small t > 0, x E [—2x /3, 2x /3] and y E (0, 1), 

\E(t,x,y)-E (t,x,y)\<Ce-* 2 o/™. 
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We observe from q < r < qe and this estimate that 

»3x /4 r l 



h,i > IT +[ }E (t,x,y)e^ t+ay (-qe y + r)dy-Ce'^ m 

\J0 J 3x /A J 

p -f3t-\a\ f3xo/4 / 1 \ 

> - / e -(y-*) 2 /*tr y + r y _ Ce -xyzu ^ = + x 

V&rt Jo J 



for small t > and x G [-2z /3, 2ac /3]. By (13TT0]) we have -ge y + r > M 2 ,i for all 
?/ G [0,3x /4] and some M 2il > 0. Thus 

h,i > M 2 for small t > 0, x G [— 2xo/3, 2xo/3] and some M 2 > 0. 

Since it is seen that ii >2 > e~ /3t_:c2//4 '/v / 8vrt for small t > and all a; G Q, we obtain 

e -/3t-x 2 /4t 

w t (t, x) > M 2 H — for small t and a; G [0, 2a; /3]. 

V87rt 

Therefore, for sufficiently small t% > 0, we have (I3.13P by integrating both sides of this 
inequality on [0,t] for all t G (0,ti). 

Step 3. Set x(t) := sup{y G [0, 1] | u(t, x) > f(x) for all x G [0,|/)} for each t G [ti,T] 

and define x 2 i := inf x(t). We claim x 2 i > 0. 

te[ti,T\ 

Suppose x 2j i = 0. Then for each n G N, there exists t n G [ti, T] such that x{t n ) < 1/n. 
Extracting a subsequence if necessary, we may assume t n — > t G [ti,T] as n — > +oo. 
Noting that u(t,x(t)) = ip(x(t)), we easily see that 

u(t, 0) = lim u(t n ,x(t n )) = lim <^(af(t n )) = <^(0) = 0. 

n— >+oo n-^+oo 

This contradicts to (I3.12p . Hence the claim of this step is proved. 
Putting x 2 := min{2x /3, x 2) i}, we obtain the desired result. □ 

Remark 3.3 It readily follows from Lemma P3 .21 that ( £ (u e — if) = in [0, T) x (— a; 2 ,a; 2 ) 
for all e > 0. Hence the boot-strap argument yields that u e G C°°((0, T) x fi). 

Based on (13. 5p and Lemma 13.21 we prove the following theorem. 

Theorem 3.4 Assume q < r < qe. Let xq be given in ( \2.3\) and x 2 G (0,xq) in Lemma 
Iff.M Then, for each x^ G (0, Xq), there exists L 3 > such that for any e > 0, 

L 3 

(3.16) ||<(t, -)\\l^(-x s ,x 3 ) < -j=, m(t,-)\\L°°(n\(-x30 a ))<L3 for all t G (0,T), 

(3-17) IImL^OIU-c-^s) < IkLC^OIU-Cnvc-^,^)) < L s for all t e (0,T). 

Proof. The u e is given by 

u £ {t,x) = u{t,x) - / / e~ l3( - t - s ' ) E(t-s,x,y)( e (u s (s,y)-p(y))dyds, 
Jo Jn 



where u is the solution of (13. lip and E is defined by ( I3.14p . Fix 23 G (0, 2 2 /2). We divide 
our consideration into three cases. 
Case 1. \x\ < x 3 . 

Differentiating u £ with respect to t, we get 



(3.18) u £ t (t,x) 



x ) - 2) - if{x)) 

e -W->)Et(t - s, x, y)Uu £ (s, y) - <f(y))dyds 



Jn 



E(t - s, 2, y)( £ (u £ (s, y) - cp(y))dyds 



'0 ./n 

" : ^2,1 + -^2,2 + h,3 + -^2,4- 



It follows from the standard theory for parabolic equations that 1/2,1! — C/^/t for all 
t G (0,T) and x G O. Besides, from we get |/ 2j2 + J 2 , 4 | < C for all t G (0, T) and 
2 G fi. As for 72,3, noting that (2 — y) 2 > x\ for all x G [—23,23] and y G [2x3, 1), we 
observe by (13. 5p and Lemma [3.21 that 



< C 



/ / E 0)t (t- s,x,y)dy + / / E 1)t (t - s,x,y)dy 

JO Jx 2 JO Jx 2 

< C I jf (t - sJ-^e-^/^-^ds + 1 J < M 3 ,i. 

for all t G (0, T), x G (— 23,23) and e > 0. Here and in the sequel, M 3ii 's (i > 1) are 
constants depending on 23. Consequently we obtain 



L°°(— £3,3:3) 



< 



M: 



3.2 



Vi 



for all t G (0, T) and small £ > 0. 



Case 2. 2 3 < \x\ < 1. 

Assume that 23 < 2 < 1. By using (13.181) and (13. 5p . it is seen that 



E t (t - s,x 3 ,y)dyds 



Jx 3 



C 



\u £ t (t,x 3 )\ < \ut(t,x 3 )\ + c 

=: 7 3jl + J 3j2 + C for all t G (0, T) and small e > 0. 
We estimate and / 3! 2- We observe by the integration by parts that 



I4il < 



E y (t,x 3 ,y)ip y (y)dy 

xs/2 



+ P\\<p\\L°°(n) 



E y (t,x 3 ,y)<p y (y)dy 
h,i,i + 73,1,2- 



x 3 /2 



E y (t,x 3 ,y)(f y (y)dy 



+ C 



Using (x 3 -y) 2 > x\/A for all y G (0, 2 3 /2], we can estimate |/3,i,i| < Ct 3/2 e < M 3j3 

for all t G (0, T) and e > 0. Since <£> is smooth in (0, 1), it follows from the integration 
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by parts and the fact \E(t,x 3 ,x 3 /2)\ + \E(t,x 3 , 1)| < Ct-^ 2 e-^ m that for all t G (0,T) 



Hence j^il < M 3 ^ for some M3 5 > 0. On the other hand, it is observed by the same 
argument as the estimate for J 2 3 in Case 1 that 



Therefore we conclude that \u\(t, x 3 )\ < M 3 7 for all t G (0, T) and e > 0. 

We provide an estimate for \\u*(t, -)\\l°°(x 3 ,i)- Differentiating the equation of (I3.1I) with 
respect to t, we have 



where U e := u\. Noting that \U £ (0, -)| < C on [s 3 , 1] and that (' £ > 0, we obtain from the 
maximum principle ||itf ||z,°°(:E3,i) < M 3>8 for all t G (0, T), small e > 0. 

We can get ||wf(t, •) II 1,-0:3) — ^3,9 f° r a ^ * e (0>^) an d small £ > by the same 
way as above. 

Consequently, for each x 3 G (0,x 2 /2), there is L 3 > max{M 3)2 , M 3j8 , M 3j9 } such that 
(gm holds for all i G (0,T). The (13371) follows from (1331) . (1331) ' and (l3X6j) . □ 

Proof of Theorem 13.21 Proposition 13.11 and Theorem 13.41 yield that there is a subse- 
quence {e n }+f^, e n \ 0, such that as n — * +00, 

v^wf", v 7 ^) — ► (Z/i, E7 2 , #3) weakly star in (L°°((0, T) x fi)) 3 . 

ByL 2 ((0,T)xft) c L 1 ((0,T)xfi), we can use (J3TTJ) to have (Ui, U 2 , U 3 ) = (u x , Viu t , Viu xx ). 
Hence the u G ^ o 7°°((0,T) x Q) follows from (J32D and these convergences. 

Set Xi G (0,x 3 ). The (13.81) and (13. 9p are derived from the above convergences and 
Theorem 13.41 The asserion of (2) follows from (13. 6p . Theorem 13.41 and Proposition 13. 11 □ 

We prepare some estimates for u xt and u tt to show Theorem 13.31 
Proposition 3.2 There exists C > such that for any small a > 0, 



and e > 





for all t G (0, T) and e > 0. 




SUp \\u xt (t, -)IU 2 (n) + |Kt||L2(( CT ,T)xfX) ^ 

te(<r,T) 



c 



(7 



Proof. Let w e be the solution of (13. ip . Then [7 £ := u\ satisfies 




in (cr, T) x Q, 
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Here the second inequality follows from (13.81) . The same argument as in the proof of [T0"1 
Lemma 3.4] yields that 

C 

sup ||< t (t, -)IU 2 (n) + IKIIl2(( ct ,t)xq) < — • 

t<E(a,T) O 

Sending e — > 0, we have the result. □ 

Proof of Theorem 13.31 Step 1. We claim that u t is continuous in (0,T) x fi. 

We observe from the regularity theory for parabolic equations that u t is continuous 
in {(t,x) | < t < T, u(t, x) > y?(x)}. It is obvious that u t is so in {(t,x) | < t < 
T, x*(t) < x < 1}. The continuity of u t in {(t,x*(t)) | < t < T} can be proved by 
Proposition 13.21 and the same argument as the proof of [9[ Corollary 4.2]. Hence we have 
the claim. 

Step 2. We show u t > in (0, T) x R 

We modify (13. ip as follows. Let {(ps}s>o be a sequence of C 2 and convex functions 
satisfying \\(pg — fWw^ia) — ¥ as 5 — > 0. We consider the following instead of (13. ip . 

f «t' £ - «i e + Pu 5 ' £ + Ce(u 5 ' £ -<Ps~ M A e) = in (0, T)_x fi, 

(3.19) < w 5 ' e (0, x) = <p s {x) for x en, 

{ u s ' £ (t, ±1) = tp s (±l) for t E (0, T) 

Here M 4 > is chosen so that (3 sup ||^||L°°(n) + C( — -^4) < 0. Then, there is a unique 

.56(0,1) 

classical solution of (I3.19P and it satisfies 

(3.20) limlim?/' £ = w locally uniformly in [0, T) x Q. 

We differentiate (13.191) with respect to t and set U s,e := u S t ' e . Then we have 

Ut' £ ~ U 5 X ' £ + (3U 5 ' £ + C(u s ' £ -<ps- M 4 e)U s > £ = in (0, T)_x fi, 
U s > £ (0, x) = tp s , xx - Pfis - C(-M 4 ) for x e f2, 

t/*' e (t,±l) = fortG(0,T). 

The [/ a ' e (0, ■) > on O follows from (ps,xx{%) — and the choice of M4. Hence we apply 
the maximum principle to obtain U s,e > in [0, T) x f2. The (13.201) and this result yield 
that u(-, x) is nondecreasing for each x £ Q. Hence we have the result. □ 



4 The discrete Morse semiflow 

Fix h > 0. As briefly mentioned in Section 2, for each m = 1,2, ... , [T/h], there is a 
unique minimizer u m G /C of the functional J m . Moreover, u m satisfies (12 .4p in the weak 
sense (cf. [3 Chapter 3]). We call the sequence {w m }m=o the DMS-BS. This section is 
devoted to some properties of the DMS-BS and to its free boundary. 

In Subsection 4.1, we discuss some properties of the DMS-BS. To prove Theorem 
14.31 of this subsection, we need an estimate of the difference u m — it m _i, Theorem 14.51 
Since its proof consists of lengthy and careful calculations, it is given in Subsection 4.3. 
In Subsection 4.2, we consider the existence and uniquess of the free boundary of the 
DMS-BS. 
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4.1 Some properties of the DMS-BS 

First, we have the monotone property of the DMS-BS. 

Theorem 4.1 Let {u m }m=d be the DMS-BS. Thenu m ^i < u m onVt for allm = 1, 2, . . . , [T/h] 
and h > 0. 

This theorem can be easily proved by the maximum principle and induction. Hence we 
omit the proof. 

In the following part of this subsection, we show the time-discrete analogues to The- 
orems 13.11 and 13.21 



Theorem 4.2 For each h > and m = 1, 2, ... , [T/h], u m is a unique solution of l \2-4 ) 
in the a.e. sense and in the viscosity sense. In addition, {u m }^[^ satisfies 



[T/h] 



sup \h ^2 



h>0 



m=l 



2 [T/h] 

+ sup |l-u m || H i(n) + h 

L 2(Q) 0<m<[T/h] m=1 



\Um,xx\\i J 2(Q') I < +00. 



Theorem 4.3 Assume q < r < qe and let x$ be defined in <\2.3ti . Then there exist 
x 4 G (0,xo) and hi > such that for each h G (0, hi), {« m }^=o C W 2,co (Q) and it 
satisfies the following estimates. 

(1) There are L 4 , L 5 > depending on rr 4 such that for each h G (0, hi) and m = 
1,2,..., [T/h], 



(4.1) 



\u m (x) - M m _i(x)| 
h 



La max 



< < 



'mh 



La 



La max 



log/i| 

u 

1 /|log/i| 



(4.2) | 

(*)l < 4 



log/i| 



^g^l l r I I 

/or \x\ < Xa, 
for \x\ > x 4 , 
for a.e. \x\ < x 4 , 
for a.e. \x\ > x 4 . 



(2) There are L 6 , L 7 > depending on Xa such that 

I M _ / xi < / L 6 (\(m - n)h\ l > 2 + \x - y\) if m, n < L 7 (h\ log h\)-\ 

\u m {x) u n {y)\ _ | L G (\(m - n)h\^f\\ogh\ + \x-y\) if m, n > L 7 (h\ \ogh\)-\ 

for all m, n — 0,1, ... , [T/h], x, y G ft and h G (0, hi). 



Remark 4.1 The -i/j \ogh\ appears in Theorem 14.31 by some technical reasons. 
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To prove Theorems 14.21 and 14.31 we introduce the penalized problem to (12. 4ft : Put 
:= if and consider 



(4.3) 



" " + /3< + ( £ (u £ m -<p) = m 0, t£(±l) = p(±l). 



Here £ £ is the same function as in Section 3. 

We observe that for each h > and m = 1,2,..., [T//i], there uniquely exists a weak 
solution of f)4.3p in the sense that u £ m — ip G #q(^) an d 



(4.4) 



' "' " '""V + "LA + W + C B « - V)0 K = for all G ^(fl). 



In addition, the regularity theory for elliptic equations yields that u £ m G H 2 (Q) n C 2 (fi). 
Thus classical solution of (14. 3p . 

We derive some uniform estimates of {Wm} m , e to prove Theorem 14.21 We get from the 
maximum principle and induction 



(4.5) 



sup i sup IKJIctn) < \m\omy 

h>0,e>0 \0<m<[T/h] ' 



By a similar argument to the proof of Lemma 13.11 and induction, we have 

Lemma 4.1 We have u £ m > ip — M\E on Q for all e > 0, m = 1, 2, . . . , [T/h] and h > 0. 
Here Mi is the same constant as in Lemma \3. 1[ 

The following estimate is a time-discrete analogue to (13. 3p . 

Proposition 4.1 We have 



[T/h] 



sup \h^2 



h>0,e>0 



m=l 



U m U m-1 



h 



L*(Sl) 0<m<[T/h] 



SUp IKJI^ 1 ^) < +°°- 



Proof. Put := in (HD- Using uu < {u 2 + u 2 )/2 = fe^-i), 

(u £ m , x ,u E m _ hx )), we have 



2 (u £ ) 2 



U 



m—l,x) 



e \2 _ ( n ,e \2 



+ (3 



m—l i 



+Ce«»-¥0«»-t4-l) 



alx < 0. 



Since it is easily seen from Lemma 14.11 that 



C h h 

Ce(u £ m ~ <p)(u £ m - U^) > -C\u £ m - > — - - 



U m U m-1 



h 
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we have 







U m U m-1 




i 5 


h 



+ 



it 



for all e > 0, m = 1, 2, . . . , [t/h] and /i > 0. Summing up these inequalities from m = 1 
to m = [t/h], we obtain 



m=l 



u e m -u e 



m—1 



h 



L 2 (n) 



+ min{l,^}||«f t/fc] ||| r i ( n) < max{l, /3}||^||| 1(n) + C. 



Since t G (0, T) is arbitrary, we have the result. □ 
From Lemma 14.11 and Proposition I4.1[ we have 



(4.6) sup 

h>0,e>0 



/ \ [T/h] 

sup \\( e (u e m -<p)\\ cm <+oo, sup /lV 

\0<m<[T//i] J h>0,e>0 m=1 



l M m,2;xlli 2 (r2) < +00. 



Proof of Theorem 14.21 From (14.5)) . Proposition 14.11 and (14. 6p . we observe that for each 
h > and m = 1, 2, . . . , [T//i], sup ||Mmll-H" 2 (n) < C/h. Hence applying Sobolev imbedding, 



£>0 



we can extract a subsequence {e n }„^, e n \ such that as n — > +oo, 



(4.7) 



7/ 



in C\Q), u 



m,xx 



U„ 



weakly in L 2 (Q), 



for all m — 1, 2, . . . , [T//i] and h > 0. Thus w m is a solution of (12 .4)) in the a.e. sense and 
in the viscosity sense. The u m = u m follows from the uniqueness of solutions of (12.41) . 
The estimates in Theorem 14.21 follows from Proposition 14.1) (14. 6 j) and (14.71) . □ 

We provide some pointwise estimates for {u £ m } mte . We can show by a similar argument 
to the proof of Proposition 13.11 that 



sup l sup \\u mx \\ L ^ {n) 

h>0,e>0 \0<m<[T/h] 



< +oo. 



The following theorem plays a crucial role to prove Theorem 14.31 and Theorem 14.71 in 
Subserction 4.2 below. 



Theorem 4.4 Assume q < r < qe and let xq be given in (\2.3§ . Then there exists hi > 
such that u £ m > <p in (—1, Xo+\^h/2) for alls G (0, h 4 ), m — 1, 2, ... , [T/h] and h G (0, h 2 ). 

To prove Theorem 14.4) we prepare some lemmas. 



Lemma 4.2 Let Mi > be given in Lemma l3A[ Then for each h > 0, m — 1, 2, . . . , [T/h] 
and e G (0, /i 4 ), we have u e m > uf n _ 1 — Mih 4 on fl. 
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This lemma is a substitute for Theorem 14. 11 Because we do not know such a monotone 
property for {u^}^ as Theorem 14.11 holds since it may happen u £ m (x) < <p{x) in view 
of Lemma 14.51 

Proof of Lemma 14.21 We can prove the case m = 1 by a similar argument to the proof 
of Lemma 13.11 Next we consider the case m = 2. Put wf := u\ — M\h A '. Since u\ is a 
classical solution of (14.31) with m = 1 and satisfies u\> Uq — M\h A on Q, we observe that 

= L T A - ul, xx + Mi + Cs(u{ -<p)< Uu\ ~ M,h A -<p)- Uu{ - y?) < in n 

and u\ (±1) < Wi(±l) for all e € (0, /t 4 ) and h > 0. Applying the maximum principle, we 
have u\ > u\ — M\h 4 on Q. . 

By induction, we obtain the desired result. □ 

Lemma 4.3 Assume q < r < qe. Then there exists hs > such that u e m > ip in (—1,0] 
andu £ m (0) > <p(0) + Vt/Afor all e > 0, m = 1, 2, . . . , [T/h] and h e (0,h 3 ). 

Proof. Step 1. We claim that u £ m > ip in (—1, 0] for all e > 0, m = 1, 2, . . . , [T/h] and 
h>0. 

For each m = 1, 2, . . . , [T//z], let U_ m be the solution of 



— m — m— 1 



+ ££m = in Q, £ m (±l) = ^(±1). 



Since is a classical subsolution of (14.31) with m — 1, we have uf > JT^ on by the 
maximum principle. We see by induction that u e m > C/ m on £7 for all e > 0, m = 
1, 2, ... , [T/h] and h > 0. Therefore, the claim of this step follows from U_ m > in Q for 
all m = 1, 2, . . . , [T//i] and (p = on [-1, 0]. 

.S^ep £. We prove that there exists > such that m^(0) > </?(0) + for all 

e > 0, m = 1, 2, . . . , [T//i] and /i G (0, /i 4 ). 

L/^ is given by 

tt ( \ 1 In ( \ ( \a j ¥>(l)sh(2 ft (a; + l)) 



h J n sh(2z h 
where := (3 + l/h, sh(r) := sinh(r) for rei, 

sh(z A (a - x))sh(z h (a + y)) 



(4.8) G aih (x,y) :-- 



(-a < ?/ < x < a) 
z h sh(2az h ) 

sh(z h (a+x))sh(z h (a-y)) 

—a < x < y < a) 



z h sh(2az h ) 

for a > and Gh '■= G\^. To estimate IZi(0), we directly calculate that for — 1 < x < 0, 

, . ~s t t / \ , , / /, hsh(zh(l + x)) ( qe a+l 

(4.9) U 1 (x) = Ash(z h (l + x))+ 1 

A 



sh(2z h ) \l + qh 1 + rh 



2/iz ft sh(2z h ) { (z h - a - l)(z h - a) (z h + a + l)(z h + a) 
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Putting x = 0, we have u\ (0) > C/^O) > Vh/3-CVhe~ Zh for all /i > 0. In view of Lemma 
14.21 and tp(0) = 0, selecting h 3 > sufficiently small, we obtain the desired estimate. □ 

The following lemma is suggested by the formal asymptotic expansion of solutions of 
(El (cf. Section 6 below). 



Lemma 4.4 Put ph := \fh — (a + l/2)h and p := (x — Xq — ph)j2\fh. We define 

f ip(x) + e ax {h 3 / 2 w 3 (p) + h 2 (w^p)- M 5 )} ioTxe[-l,x + ti h ], 
UX> ' \ <p(x) - M 5 h 2 e ax for x e (x + p h , 1] 

(4.10) w 3 (p):=r(e 2p -l-2p), 

w 4 (p) := r{e 2f > - (1 + 2p + 2p 2 ) + a(e 2p - 2pe 2p - 1)}. 

Then there are large M 5 > and small > suc/i u is a subsolution of with 
m = 1 in the a.e. sense satisfying w(±l) < </?(±l) /or a// £ G (0, /i 4 ) and h G (0, /14). 

Proof. Note that u G W /2 '°°(fi) fl C 2 (Q\{x + /i h ,}) and that w 3 , u> 4 satisfy 



(4.11) 
(4.12) 



w 3 - ~f + r{2p + 1) = 0, W4-y-a^ + r{2p 2 + 2p - a) 
\Vhw 3 \ + \hw 4 \ + \Vhw' 4 \ <C in (-(1 + x + p h )/2Vh,0}. 



We divide our consideration into two cases. 
Case 1. x G (— 1, £0 + jU/J. 

In this case, p G ( — (1 + xq + ph)/2vh,0]- Using (14. lip , we compute that 



U — (f 

h 



u, 



+ (3u + ( £ (u -tp) < e ax [h{-M 5 + r(Vhw 3 + hw 4 ) - aVhw' 4 } 

+qe x -r - Vhr(2p + 1) - hr(2p 2 + 2p - a)}. 



We see from x = log(r/g) and x — xq = pu + 2\/hp that 



qe x — r 



r e r 



(// h + 2Vhp) + ^(p h + 2Vhp) 2 + h 3 ' 2 



< r{Vh(2p + 1) + h(2p 2 + 2p-a)} 

1 - Vh(2p + 



+rh 



1 

a+ 2 



Here we have used ph + 2\fhp < \/h/2 and the following inequality: 



3* - 1 - f 



if f < 0, 
h 3 ' 2 if £ G [0, 



for small h > 0. 



By the fact -y/h(2p + 1) < Xq, we get 



- r - r{Vh(2p + 1) + /ir(2p 2 + 2p - a)} < C7i. 
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From f)4.12p and this estimate, we have 



u — (p 



u xx +/3v + Ce{u -<p)< he ax (-M 5 + C) in (-1, x + p h ] 



h 



for small e > and h > 0. Taking M 5 > large enough, we conclude that u is a classical 
subsolution of (14.31) with m — 1 in (—1, x + Ph) for small e > and h > 0. 

Case 2. x & (x + ph, 1). 

Taking M5 > 1 and small h$ > 0, we see that for all e G (0, /i 4 ) and ft G (0, ft.5), 



Therefore for large M5 > and small > 0, tt is a subsolution of (14. 3 p in the a.e. sense 
for e G (0, /i 4 ) and h G (0, /14). In view of (14.121) . we can get u(±l) < y(±l) by replacing 
h± with a smaller one if necessary. Thus the proof is completed. □ 

Proof of Theorem 14.41 It follows from Lemma 14.41 and the maximum principle that 
u\ > u on Q for all e G (0, h 4 ) and h G (0, ^5). In view of Lemma [4.31 we have only to 
prove the assertion on [0,Xq + y/h/2]. 

First we treat the case m = 1. Let p and fih be defined in Lemma 14.41 and set 
Pi := — (xq + Hf^/^Vh, P2 '■= —1/4 + (2a + l)y/h/A. We observe by careful calculations 
that for small h > 0, 



Hence we have h 3 ^ 2 ws + h 2 (w4, — M 5 ) > r/i 3 / 2 /20 on [pi, P2} and thus 

u\ > <p + —h 3/2 on [0, xq + Vh/2] for any e G (0, h A ) and small h > 0. 
20 

In the case m > 2, Lemma [4.21 and the above estimate yield that 

u e m > up + -U?l 2 - MiTh 3 > V + ^h 3 ' 2 on [0, x + Vh] 

for all e G (0, h A ), m = 1, 2, . . . , [T//i] and small h > 0. 

Hence, selecting h 2 > sufficiently small, we have u e m > <p + rh 3 / 2 /40 on [0,a;o + vh/2] 
for all e G (0, /i 4 ), m — 1,2, ... , [T/h] and ft, G (0, h 2 ). Thus we complete the proof. □ 

By Theorem I4.4[ we see that for any h G (0, hi), m = 1,2, ... , [T/h] and e G (0, h A ), 
u e m satisfies u e m > <p in (— x , x ) and thus 



u — ip 



h 




X 




J^{h 3/2 w 3 + h 2 (w 4 - M 5 )} < on [ Pl , p 2 ], 
h 3/2 w 3 ( Pl ) + h 2 (w,( Pl ) - M 5 ) >^h- M 5 h 2 , 
h 3/2 w 3 (p 2 ) + h 2 (w 4 (p 2 ) - M 5 ) > ^-h 3 ' 2 - M 5 h 2 . 




h 



m,xx 



+ (3u e n 



rn 



in (-x ,x ). 
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We prove Theorem 14. 3^ based on this fact. Before doing so, we give some preliminary 
analysis. 

In (-xq,x q ), u e m is given by 

<(*) = \f G n ^y)u' m _Mdv+ u " ( ~ Xa)sHzh(x, '~ x)) ■ 



h in ° m ~ ly ' sh(2x z h ) sh(2x a z b ) 

where G XQi h is defined by (14.81) with a = Xq. In the sequel we set Xq — 1 for simplicity. 
Define ^ 

W]{x) :=lj G h (x,y)i/>(y)dy for ^ e C([-l,l]). 

For — 1 < x < 0, we get from (14.91) 

hsh(z h (l + x)) ( qe a+1 re a \ y9(l) sh(z fe (l + a;)) 

SfcM(x) = sh (2 ,,) (tt^ ~ ttW sm^j — + AshM1 + 

On the other hand, we observe by tedious calculations that for < x < 1, 

/ish(3 A (ac + 1)) f ge a+1 re c 



sh(2z fc ) \l + g/i 1 + r/i 



2hz h sh(2z h ) { (z h + a + l)(z h + a) (z h — a- l)(z h - a) 
Noting that (p = on [— 1, 0] and that A > B > 0, we have the following: 

9 h [<p){x) < <p(x) + R(x) - ^ 1)S ^ h) +X)) for all x G [-1, 1], 

_ hsh(z h (l +x)) f qe a+1 re a \ / ' qhe {a+1)x rhe ax \ 

(x) : ~ sh(2^) [TT^h ~ TTVK ) ~ { i + qh ~ TTT^ ) 1{x>0} 

+Ash(z h (l-\x\)). 

Here l{ x> oy(x) = 1 for x > and = for x < 0. Recalling Uq = (p, we see that 

ul{-l)sh(z h (l-x)) , uf(l)sh(z h (l + x)) 



sh(2z fc ) sh(2z fc ) 



< u (x) + R(x) + Mf( 'l/f 1} sh(* fc (l - *)) + " !(1 ^ 9 "° (1) sh(„(l + x)). 
We can inductively show that 

<(x) < <,-iM + sr'ww + E "' t( ~ 1 ip^)' ( ~ 1) srt[sh(z ' i(1 ~ 
+S%^ ,( V'i*fa(i+-))i(»), 

where S|[^] := S^fS* -1 ^]] and 9/1I/0] := V ; - Thus we need some pointwise estimates for 
5Z,hi R ] and $2,h[ sh M x ° ± •))] to prove Theorem O 
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Theorem 4.5 For h > 0, m — 1, 2, . . . , [T//YJ and x G [— x ,x ]; we ^aue 



(4-13) 9™^[sh(%(xo - I ■ |))](s) < sh(z,(x - |x|)) ^(^\x\) k + Ch 



k=0 



k\ 



m , 

(4-14) 9SU[ah(s fc (so ± •))](*) < Bh(^(x ± *)) £ — t^(^(x T *))' 



2m — k k\ 

fc=i 



w/iere a mifc := (2m — fc)!/{2 2m k m\(m — fc)!}. 



Theorem 4.6 There exist x^ G (0,Xo) and L 7 , L 8 , L 9 > and /i 5 > such that for each 
h G [0,h s ), 

(4.15) ^ h [R]{x) < { forxe(-x 4 ,x A ), fora ii m = h2 ,...,[T/h], 

y L 7 h for x = ±x 4 , 



,4,6) ■))](*) 

/or aZZ m = 1, 2, . . . , [Lg/h\ logh\] and x G [—2:4, X4]. 



We admit that Theorem 14.51 holds and prove Theorems 14.61 and 14.31 We give the proof of 
Theorem 14.51 in Subsection 4.3 below. 

Proof of Theorem 14.61 Set xq = 1 for notational simplicity. Note that 

(J e -k 2 /4(2m-k) 

(4.17) a m)k < — -p===— for A; = 0, 1, ... ,m. 

This will be proved in Subsection 6.2 below. 

Step 1. We show (I4.15P for some x^i G (0, 1) and L 7 > 0. 
Since it is easily seen from (14. 13j) that 

m 

SJM(aO < V^c"* sh(z fc (l - |x|)) ]T 2=h*(* fc |s|)* + C%, 

fc=0 



we have only to treat the first term of RHS of this inequality. Denote it by I±{x). 

h\x\ 
k\ 



m ( I \\k 

Using (147171) and ^ !, < e* fe|a:| , we get 



fc=0 

(4.18) J 4 (x) < CV— for all m = 1,2,..., [T//i], xG [-1,1] and /i > 0. 

V 771 

Fix x^\ G (0, 1). We consider only 14(2:4,1) since I4 is even. We still denote it by J4 if 
no confusion arises. Set mi := [2zfc,X4,i/5] and m 2 := 3[^X4,i]. Then we have from (14.171) 

m ~fc 2 /4(2m-fc) 1 
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We divide our consideration into three cases. Let ft 2 be given in Theorem 14.41 
Case 1. m < m\. 
Using Stirling's formula, we get 

Set 7 = k/z h x 4 ,i. Then 7 G (1/^X4,1,2/5) and (z h x 4 ,ie/k) k = exp(^x 4 ,i(-7 log7 + 
Since we see that —7 log7 + 7 < 4/5 for all 7 G (1/2:^4,1, 2/5), we have 

Hence J 4 < Ch^e'^ 1 / 5 for all m = 1, 2, . . . , [T/ft] and h G (0, ft 2 ). 
Case 2. m> m\. 

We may assume m > m 2 . Similar calculations as in Case 1 yield that 



mi 7712 



h < CVhe~ ZhX ^ [y^ + + Yl 



-fc 2 /4(2m-fc) 1 

1 \fc 



K fc=0 fc=OTl+l fc=7Tl2+l/ 



V2m-k k\ 



' ' ' - e -fc 2 /4(2m-fc) j 



e -fc 2 /4(2m-fc) 

K=7711 + l 

Setting m := l/ft,s, we observe that for k — mi + 1, ... , m2, 

-fc 2 /4(2m-fc) /m 

- < e - m ?W8 < cVhse~ M4 ' lS < M^\[h. 

\/2m-k s/m 

Here and in the sequel M 4)i 's (i > 1) are positive constants depending on x 4j i. Hence we 
have 

I A2 <M A2 Vh V < ^4 2 Vhe ZhX4 - 1 for all ft G (0, ft 2 ). 

z — ' re! 

fc=771l+l 

Consequently we obtain 7 4 < M 4i3 ft for all m > mi and ft, G (0, ft 2 ). 

From Case 1 and 2, choosing L 7 large enough, we get (14. 15j) for all ft, G (0, ft 2 ). 

5"tep 5. We show that for any X4 i2 G (0, 3/4), there are h^ t \ > 0, L 8 , Lq > such that 
fT4~T6]) holds for all m = 1, 2, . . . , [L 9 /ft| logft|], X G [ — X4 j2 , X4 ;2 ] and ft G (0, ft.5,1). Denote 
by h{x) the left-hand side of fl4T6l) . 

Fix x 4j2 G (0,3/4). For x G [— x 4j2 , x 4j2 ], set m\ := [2^(1 — x)/5], m 2 := 3[^(1 — a;)] 
and m 3 := [^(1 - x) 2 /100| logft|]. We use (PI]) and sh(z h (l + z))/sh(2z A ) < 2e- Zh{1 ^ 



on [—1, 1] to have 

rn 

h(x) < Ce- Zh{l ~ x) Y^ 



g Ze -/V4{2(m-fe)-0 ( Zfc (l_a;))< 



fe=l 

m I m 



f {2(m - k) - iy/ 2 

, e -t*/4(2fe-0 I ( ^ ( 

(2fc - Z) 3 / 2 I /! 



=1 k fc=/ 
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We divide our considerations into two cases. 
Case 1. m < mi. 

It is easily observe from the fact le~ l2 /A{ - 2k ~ l) / \?2k - I < C for all k, I G N that 

m ] e -l 2 /H2k-l) m i rT/h i 

Epfc^^E^^ 5*^"**l foralW = l,2,...,m. 

Hence we use this inequality and the same argument as in Case 1 of Step 1 to obtain 

h{x) < C\ \ogh\e- Zh{1 - x) J2 ^-^y < C \ \ gh\e- Zhi - x - x)/ s for all x G [-z 4 ,2, x 4)2 ]. 

i=i 

Case 2. mi < m < m 3 . 

We may consider m > m 2 . Similar calculations to Case 1 yield that 

<- ^-*> E+ E + E E ';,W J)) 

\/=l i=mi+l Z=m 2 +1/ fc=Z V 7 

J" 2 , ™ ]p -l 2 /4(2k-l) l 7 H _ T \\l 

< dJogfcle^-W+Oe-^-) E E L _ A3/2 ,, — ■ 

m \ i c ~i 2 /i(2k~i) / m x \y 

Set J 5 i0) = V > 7771 ,xo /0 7i — — • From the facts Z < 2k - I < 2m 3 and 

i= mi +i k=i y ' 
I > mi, we see that for k = I, . . . ,m and I = mi + 1, . . . , m 2 , 



/ — l 2 /4(2fc— /) i 

l ±- < ^ e ~-?/8-3 < C/l l/4 e -2|logh| = ^9/4 

(2/c — I) I 

Thus for x G [— 2 4j2 , £ 4 , 2 ] and small /i > 0, 

/ 5 ,i(x) < c jr jrh 9 / Azh{1 - x))l <ch^m 3 jr MI_e02! 

l=m,\-\-\ k=l i=mi+l 

< M 5jl /i 5/4 | log/ile 2 ^ 1 "^ < M 5A he Zh(1 - x) , 

where M 51 depends on 1 — x 4)2 . Consequently, we get 

h,i( x ) — M 5i ih for all x G [— x 4)2 ,x 4)2 ] and small h > 0. 

Thus taking large L$ > 0, L 9 := (1 — x 4j2 ) 2 /100 and /15 := min{/i 2 , /i 5j i}, we obtain 
flUED for all m = 1, 2, . . . , [T//i], x G [-x 4 ,x 4 ] and /i G (0, h 5 ). 
Setting x 4 := min{a; 4jl , x 4j2 }, we complete the proof. □ 

Proof of Theorem 14.31 Step 1. We claim that there are L 4 ,i, £5,1 > and hi > such 
that for all h G (0, hi), m = 1, 2, . . . , [L 5i i//i| log and x G x 4 , x 4 ). 

(4.19) u E m {x) - u^x) < L 4jl h. 
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Put hi := h e . It follows from Theorem 14.61 that 



(4.20) <(±X4) - <-l(±^4) < (^7 + ^ 8 )/» 

for all m — 1, 2, . . . , [L g /h\ logh\] and h G (0, foi). Choosing L 4 > L 7 + L 8 , we observe 
that Uq + L 4 /i is a classical supersolution of (14.31) in fi\(— £4,0:4) since Mq(= V 9 ) is smooth 
in this domain. Hence we apply the maximum principle to have (14.191) with m = 1. We 
inductively obtain (l4.19p for m = 2, . . . , [Lg//i| log/i|] and h G (0, /ii). Putting L 5 := L 9 , 
we have the claim. 

Step 2. We derive the estimates of (1) and (2). 

From Theorem 14.61 and (14.191) . we obtain 

if w < .„ 

V m 

L4/1 if b| > x 4 



(4.21) u^r)-^^) < 

for all m = 1, 2, . . . , [L5//1I log/i|] and small h G (0, hi). 



Hence OH]) holds for all m = 1, 2, . . . , [L 5 //i| log/i|]. 
In the case m = [L 5 /h\ log/t|], we have 



lOg I 

«m - <-l < ^4^W ' Ls On fl. 

We can show by the maximum principle and induction that this estimate holds for m = 
[L 5 /h\ logh\] + 1, . . . , [T/h]. Using (14. 7p and the above estimates, we have (14.11) for all 
m = 1, 2, . . . , [T/h] and h G (0, h x ). 

The (14. 2 p can be derived from (14. ip . (14. 5p and (14. 6p . The estimate of (2) is a conse- 
quence of IjOl and (JUTJ). □ 



4.2 Free boundary for the DMS-BS 

The problem (12.21) has a unique free boundary. However, it does not leads to the existence 
and uniqueness of that for the DMS-BS. To prove them is the purpose of this subsection. 

Theorem 4.7 Assume q < r < qe. Then there is h§ > satisfying the following: For 
each h G (0, h 6 ) and m = 1, 2, . . . , [T/h], there exists a unique x m G (xq, 1] such that 

(4.22) {x eVL\ u m (x) > if{x)} = (-1, x m ), {x G | u m (x) = (p(x)} = [x m , 1]. 

Moreover, Xq < Xi < x 2 < ■ ■ ■ < x m < ■ ■ ■ < x\T/h] an d x m < min{xo + \fmh, 1} for all 
m = 1,2,..., [T/h]. 

Proof. Put h 6 := minj/i!, h 2 , h 5 }. Notice by Theorems 14.11 and 14.41 Lemma I4T31 and (14.71) 

that u m > (p in (— 1,xq + Vh/2) for all h G (0, h§). Since we easily observe by Theorem 
14.11 that x\ < X2 < • • • < x m < • • • if they exist, in the following we show the existence 
and uniqueness of x m . 
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Step 1. We treat the case m — 1. Set p = (x — x — V~h)/2Vh and 



uAx) :-- 



(f{x) + e ax h 3 / 2 w 3 (p) for x G (0, x + v^), 
(f(x) for x G [xq + Vh, 1), 



where w 3 is defined by (14.101) . We show that u\ is a supersolution of (12 .4p in (0,1) in the 
a.e. sense. 

Note that u x G W 2 <°°{tt) fl C 2 (£l\{x + Vh}) and that ¥i > <p on [0, 1]. We see by 
flUIID that in (0, x + Vh), 

- ui, xx + put = e ax {-haw' 3 + rh 3/2 w 3 + r(e^ (2p+1) - 1 - Vh(2p + 1))}. 
It follows from the facts w 3 , —w' 3 > in (— oo, 0] and e y > 1 + y for all y G R that 



h 



- u 1)XX + (3ui > in (0, x + Vh). 



By (13.101) and this inequality, we see that ui is a supersolution of (12.41) in (0,1) in the 
a.e. sense. 

In view of ui(0) < ui(0), we modify u\ to construct a viscosity supersolution of (12.41) . 
Put r] := 2(x + Vrf. Define W l by 



Wi(x) 



-7(x - xofe-^l^-^ 2 (0 < x < x ), 
(x < x < 1), 



where 7 >_0 is selected_later. Then W x G C 2 (0, 1), > 0, W liX < in (0,1) and 
W\(xo) = Wi tX (xo) = Wi tXX (xo) = 0. Moreover, it is easily observed by the choice of 77 
that 

(4.23) - W ltXX + 0Wi > in (0,1). 

Take 71 > satisfying W\ > ||y?||L°c(n) on [0, xq/4]. Since W 2 '■= |MU°°(n) is a classical 
supersolution of (12.41) in Q, setting 



UAx) 



W 2 if - 1 < x < 0, 

min{wi(x) + Wi(x), W 2 } if < x < 1, 



we conclude that U\ is a viscosity supersolution of (12.41) satisfying f/i(±l) > <^(±1). 

We have «i < U\ on [0, 1] from the comparison principle for viscosity solutions. Using 
this inequality, we can obtain a unique X\ satisfying (14.221) . Indeed, U\ = (f on [xo + V^, 1] 
implies that ui = (f on [x + Vh, 1]. Put 

X\ := mi{y G f2 | U\{x) = <p(x) for all x G [?/, 1)}(< Xq + Vh). 

Clearly u\ = ip on [xi, 1]. To verify u\ > (p in (xo + Vh/2,x\), we suppose that there is 
xi G (xo + Vh/2,xi) such that «i(x) > y?(x) in (x"i,Xi) and Ui(xi) = <£>(xi). Since U\ and 
(p are solutions of (12 .4p with m — 1, we get «i = on [x"i,xi] by the uniqueness. This 
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contradicts to the definition of X\ and hence U\ > (p in (xo + \fh/2,Xi). This observation 
also leads to the uniqueness of X\. Therefore we have the desired result of Step 1. 
Step 2. We prove the case m = 2. 

Let W3 be defined by (I4.10p and set p m = (x—XQ — \/mh)/2\/mh for m = 1, 2, . . . , [T/h]. 
By the facts w' 3 < on (—00, 0] and p m < p m _i < on [0, xo + y/(m — l)h], we see that 

(4.24) w 3 (p m ) > iu 3 (/W-i) for all x G [-l,x + \fijn - l)h] and m = 1,2, . . . , [T/h]. 
Define 



w 2 (x) := 



(p(x) + e ax {2hfl 2 w z (p 2 ) for x E [0, x + V2h], 
fix) for x G [xq + V2h, 1). 

W 2 i-l<x< 0), 

min{w 2 (x) +TVi(x), TF 2 } (0 < x < 1), 



We claim that U 2 is a viscosity supersolution of (I2.4p with m = 2. 
It is observed by (14.231) and Mi < tti + W\ in (0, xo + vh] that 

r C/a,** + P^2 > 7 v>2,xx + pu 2 

h h 

> e ax [Vh{(2V2 - l)w 3 (p 2 ) - wM} - 2haw' 3 (p 2 ) 
+r(2h) 3 / 2 w 3 (p 2 ) + r { e ^(2p 2 +i) _ 1 _ V2h(2p 2 + 1)} 

in (0, Xo + Vh}. Using P~TT]h ffl~24D with m = 2 and > 1 + y for all y G R, we get 

^ 2 ~ Ml - F 2 ,^ + ff^a > on(0,x + v^]. 

From f/ 2 > <p on ^) (I3.10P and this inequality it follows that U 2 is a viscosity supersolution 
of ([23D with m = 2. 

Thus we use the comparison principle for viscosity solutions to obtain u 2 < U 2 on Q. 
We can show by a similar argument to Step 1 that there exists a unique x 2 satisfying 
( Q2]) with m = 2 and x 2 < x + V^. 

Step 3. We consider the case m > 3. 

By induction we assume that there exists a unique x m _i satisfying (I4.22p and x m -\ < 
xq + \J {m — l)h. Then define 



Um\X 



Um(x) 



(p(x) + e ax (mh) 3 ^ 2 w 3 (p m ) for x G (0, x + y/mh], 
<f(x) for x G [xo + Vmh, 1). 

W 2 (-1 < x < 0), 

min{?I m (x) + Wi, TF 2 } (0 < x < 1). 



By a similar argument to Step 2, we can see that U m is a viscosity supersolution of (12.41) 
and thus u m < U m on Q. Therefore by the same way as in Step 1 we can find a unique 
x m satisfying (I4.22p and x m < min{xo + Vmh, 1}. □ 
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4.3 Proof of Theorem 14.5 



To prove Theorem 14.5} we prepare some identities. Define 

(z h y) k ch(2z h (l+y))dy, 3 k ■ = / (z h y) k ch(2 z h y)dy, 



i 



X k := J (z h y) k di(2z h (l-y))dy, H k := J (z h (l - y)) k ch(2z h (l + y))dy, 
M k := [ (z h (l - y)) k di(2z h y)dy (k = 1, 2, . . . , m). 



Direct calculations yield that 

D = — gk(2z h ), \ = ^—(ch(2z h ) - 1), 

do = ^-sh(2z h x), = 2~ |A,ish(2^x) - P 2}1 ch{2z h x) + X - 

%o = ^- sh(2^(l - X )),X 1 = -L {P 1X sh(2^(l - x)) + P 2)1 ch(2^(l - x)) - P ,i} , 
2^/i 2^ 

£ = 7T- sh(2z ft (l + a;)), ^ = J- {Qi.i sh(2^(l + x)) + Q 2)1 ch(2* h (l + x)) - Q 0) i} , 
2^ 2^ 

M = -^(sh(2z h ) - sh(2^a;)), M x = ^ {Q 1)1 sh(22 fc a;) + Q 2)1 ch(2z^) - Q 0jl ch(2z h )} 
2z h 2z h 

where P ,i := 1/2, P M := ^x, P a ,i := 1/2, Qo,i := 1/2, Qi,i := z ft (l - x), Q 2 ,i := 1/2. 
For k >2, the following identities hold. 

Lemma 4.5 Let J^, 3k, ^k, &k, fee defined as above. For k > 2, we have 
(-1)% = ^ (-iV + %^sh(2„) + ^ch(2z fc ) 

dh = 2^~ h ( P i^ sh ( 2 ^) - iVch(2z h x) + , 
k} 

% k = —(P ltk sh(2z h (l - x)) + P a>fc ch(2z fc (l - x)) - P ,fe), 
2z ft 

fcl 

£fc = — (Qi, fc sh(2^(l + a;)) + Q 2 , fe ch(2^(l + x)) - Q , k ), 
M fe = — |(5i ifc sh(2^x) + Q 2jk ch(2z h x) - ^^-sh(2z h ) - ^0^-ch(2z h ) 



where 



0,fe ' ^2 l (k-l)V 1,fc ' 2^2 l (k-l)V 2 ' k ' 2^2\k-l)V 

I :odd i :cvcn Z:odd 

n >f (2^) fc -' n A (**(!-*))*-' n _ A(z fc (l-x))^ 



tS'^Jfe-O! ' 2 l (k-l)\ ' w ' ^ 2 l (k-l)\ 



1=0 

Z:even Z:odd 
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Proof. Integrating by parts we have, for k > 2, 

(-1)% = - \ kz i~ l + - i)(-i) fe ^-2 ^ 

a fc = |(^x) fc sh(2z fc x) - ifc(z h x)*- 1 ch(2z h x) + - l)J fc _ 2 

X fc = j (^x) fc sh(2^(l - x)) + ^(^) fc - 1 ch(2^(l - x)) - ^(-za)*" 1 

+ ^(A; - l)£, fc _ 2 

M h = -L j - (z fc (l - x)) fc sh(2^x) - ~(s fc (l - x)) fe - 1 ch(2^x) + - l)M fc _, 

Using these recurrence formulae, we obtain the result. □ 

We separately prove (14. 13j) and (14.14p of Theorem 14.51 Put Xq — 1 for the sake of 
simplicity. 

Proof of (gZEfl . Set G m (z) := S^[sh(^(l - | ■ |))](x) for x E [-1,1], sh(r) := sinh(r) 
and ch(r) := cosh(r) for r6l Note that 



hz 2 h 



f 1 1 

(4.25) G h >0 in (-1,1) x (-1,1), J G h (x,y)dy < 

In this proof we use the identities in Lemma 14.51 and the following ones. 

(4.26) sh 2 (^(l + y)) = ~(ch(2^(l +y)-l), 

(4.27) sh(2 h (l - y))sh(z h (l + y)) = \{ch{2z h ) - ch(2z h y)), 

(4.28) -sh(z h (l-x))sh(2z h x) + sh(z h (l + x))sh(2z h (l - x)) = sh(2z h ) sh(z h (l - x)), 

(4.29) sh(* fc (l - x))ch(2z h x) + sh(z ft (l + x))ch(2^(l - x)) = sh(2^)ch(^(l - x)). 

Step 1. We estimate G\(x). 

We calculate with using (jOBl) and P~2"7|) to get for < x < 1, 

+*)) l.^^.^. 



2/iz/ l sh(2z; l ) 2z ft 
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By g2BD, we get 



where </i(:r) := (1 + z^x) sh(^(l — s)) and b := ch(2z^)/ sh(2z^). By similar calculations 
we have G\(x) < bgx(—x)/2hzf l for x G [—1,0]. Hence we obtain 



(4.30) 



Gt(x) < — ^(Izj) for 16 [-1,1]. 



U. We consider the case m — 2. 
It directly follows from (1425]) and (OTO that G 2 (x) < bS h [gi)(x)/2hz 2 h for a; G [-1, 1]. 
We observe by (gj26D, (B~2T|) and Lemma S3] that for < a; < 1, 



Sh\gi](x) 



sh(z h (l - x)) 
2hz h sh(2z h ) 

, sh(^(l + x)) 



2/, 



3 + (-l)Ji + — U fc ac + irXzhxy } ch(2z h ) - {do + 3i} 



2! 



2hz h sh(2z h ) 
-■■ /:,:• 



We get from g2HD and (14291 



4i 



2hz 2 h sh(2z h ) 



-{2 + z fc a;}sh(z fc (l - x))ak(2z h ) 



1 



+ ^ + z h x + ^(^x) 2 [> sh(z fc (l - x))(2ch(2^)) 



+- {ch(z h (l - x)) sh(2z h ) - sh(z ft (l + x))} 



Since direct calculations yield that for all x G [0,1], 

(4.31) ch(^(l - x)) sh{2z h ) - sh(^(l + x)) < sh(z h (l - x)) sh(2z fc ) + -e~ z \ 



we conclude that 



5h[gi](x) < ^92(x) + T 2 for < x < 1, 



where 



02 (a) := 



sh(^(l — x)) (3 3 



2/^ 



+ -z h x + — (z fc x) ^ , 7 



2 2 



2! 



8(2hz 2 h )sh(2z h y 



We can obtain by the similar way as above 9/i[pi](a;) < bg 2 (—x)/2hz 2 ,+'?2 for — 1 < x < 0. 
Consequently, we have 



G 2 (x) < 



(2hz 2 ) 2 



g2{\x\) + b?2 for xe [-1,1]. 
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Step 3. We estimate G m (x) by induction. 
We assume that 



b< 



m—l 



m—l 



G m -i(x) < . g m -i(\x\) + J2 h ^ fOT x e !]» 



(2^ 

where <7 m _i and 3^ are defined by 

/ m— 1 

g m -x(x) := Yl 



1=2 



( '"- l '\z h x) 1 I sh^(l-x)), '.)>/ : 



l-l 



1=0 



/! 



4(2^)'sh(2z fc )^ 



t 2q (P 

g:odd 



We estimate 5h[gm-i](x) for < x < 1. We calculate that 



2^sh(2^) ^ fc! v 

sh(z ft (l + x))^c m _i ifc 
+— rr^ — — — n — A fc 



ch(2z h ) - a* 



2hz h sh(2z h ) f^ Q k\ 



+ ^2D l ~ 1 T l forO<i<l. 



1=2 



It follows from Lemma [4.51 that 

J 5 , 2 := sh(z h (l-x)){(-l) k J k -dk}+ sh(z h (l+x))X k 



k\ 

2z~ h 



sh(z h (l - x)) { -P , k + ^Psh(2z h ) + ^ch(2z,) 



-P ljk sh(2z h x) + P 2tk ch(2z h x) 



Lfc:odd 
2 k 



+ Sh(*fc(l + x)) {Pi, fc sh(2z ft (l - x)) + P 2)k ch{2z h {l - X)) - P , k } 



Using (TOE]), (g^HD and sh(2* h ) < ch(2z h ), we obtain 
fc! 



h,2 < 



2z h 



sh(z h (l - x)) |p llfc sh(2z fc ) + ^ch(2^)| + P 2 ,feCh(^(l - z))sh(2;z h ) 



-sh(^(l + x))P ,, 



From f)4.3ip and the fact P 2 ,fe < Po,fc on [0, 1] we get 
k\ 

h,2 < 



2z h 

Consequently, we have 



sh(^(l - x)) { (P 1>k + P 2tk )sh(2z h ) + ^^(2^) } + e- z »P 0ik 



bsh(z h (l - x)) 



m—l 



2hz\ 



-Zh 



k=0 
m—l 



5/,[//,„-i](.r) < — 1 C ™ u ' {{[^ ' ! 



fc+i 



2 l {k-l)\ J (Jfe + 1)! 



fe-i m -! 



4.2^sh(2z h ) 



fc=0 
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Therefore setting 



(4.32)<7 m (x) := sh(z,(l - x)) ]T 

fc=0 

= shW i-x))x;^U(^ + ^^7v) + 



fc=0 

we obtain 



i ( i , 

2-<y{k-i)\) + (k + iy. 



G m (x) < -— ^— (/ ro (z) + ^fr*" 1 ?* for x G [0, 1]. 
K znz h) k=2 

b m m 

Since we see G m (x) < 2 , g m (-x) + yjZ/ 1-1 ^ for — 1 < x < by the same way 

^ ^ fc=2 

as above, we conclude that 

G m {x) < p^UN) + X> fe - 1 ? fe for x G [-1,1]. 



fe=2 



.S'tep ^. We determine {c m) fc}£L for to = 1, 2, ... , [T//i]. 

From (14.321) . we can obtain the following recurrence formulae: for m — 2,3, ... , [T/h] 
and k — 2, 3, . . . , to — 1, 



_ 1 

fc to— 1 

Cm— 1,1 



(4.33) Cm^jn—k — ^ ^ 9A;— z ^m—l^—l—l J Cm,o — ^ ^ 

First, we easily get 



z=o z=o 



c m , m = 1, c m _ m _i = ~(to + 1) for to = 1, 2, ... , [T/h]. 
As for c mim _2, using (14.331) and these formulae, we have 

c m , m - 2 = 7p-y( m + l)[m + 2). 

We assume by induction that for m > 3, 

i 

1 n 



^yJJ(m-l+p) for Z = 1,2,..., to- 1. 
' p=i 

From (I4.33P and this equality we compute that for = 2, 3, . . . , to — 1, 

Z=0 ' p=l I ' p=l 1=3 ' p=l ) 

lfl 3 - 1 ' 1 

l ' p=l Z=4 ' p=l J ' p=l 
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Cm,m—k 



Consequently, replacing k with m — k, we obtain 

(2m - k)\ r , 

c m fc = ; — T7 ttt for k = 0, 1, 2 . . . , m. 

' 2 m - k m\(m - k)\ ' ' 

Step 5. We derive fQ5]) . 

It is easy to see that for small h > and m = 1, 2, . . . , [T//t], A < (1 + Ce~' iZh )/z h e Zh 
and (b/hzh) m < 1 + e~ 22h . Besides, since c mj fc/2 m is the m-th term of the binomial 
expansion of (1/2 + l/2) 2m_fc , it is obvious that c m ^/2 m < 1. Using these facts, we get 

snsh^(i-i-D)](x) < — um) + -^E|+^+E^ 



< 

z h e Zh 



^ m 

g m (\ X \)+Ch + J2^l 



for small h > 0. Similarly we observe that 



m ^-y _ 2 m i-1 p / \„ m Z—l -, 

E^^EE^E^^ EE ^<c A -v 

i=2 v y Z=2 p=0 9=0 * Z=2 p=0 



for all m = 1, 2, ... , [T//i] and /i > 0. Setting a mj fc = c mj fc/2 m , we obtain (I4.13p . □ 

Proof of ( 14.141) . We treat only 9™[sh(zk(H-))] because 9/j l [sh(2:/ l (l — ))] can be similarly 
estimated. Set H m (x) := S^[sh(^(l + -))}(x). In this proof, we use the indentities in 
Lemma 14. 5 [ (14.261) . (I4.27P and the following ones. 

(4.34) sh(z A (l -a0)sh(2z h (l + a;)) - sh(z h (l + x))(sh(2^) - sh(2z fc x)) = 0, 

(4.35) sh(z fc (l - z))(ch(2z h (l + x)) - 1) + sh(^(l + x))(ch(2z h x) - <h(2z h )) = 0. 

Step 1. We consider the case m — 1. 
Using (14.261) and (I4.27p . we compute that 

. , sh(zfe(l — a;)) 1 , , 



From (I4.34p we have 

(4.36) #i(z) < —3/11(2;) for are [-1,1], := - x)sh(z fc (l + x)). 



Step 2. We estimate the case m = 2. 
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It follows from (TOoD and (TOoD that H 2 (x) < b9 h [h 1 ](x)/2hz 2 h for -1 < x < 1. We 
see from P~2T)j) . fTOTl) and Lemma S3] that 



Sfc[/ii](a:) < 



sh(z /t (l - x)) 
2/iz /l sh(2z /l ) 

| sh(z fe (l + x)) 
2fe ft sh(22 /l ) 



1 

2^ 
1 



(z h (l - x)) sh(2^(l + x)) + -(ch(2z ft (l + x)) - 1) 



2 j2fc(l - ^) 2 ch(2^) + 2fc(l - x) sh(2^a;) 



1 



+ -(ch(2^x) -ch(2z h )) 



We use fljjjgp and (1435]) to obtain 



2/^ 



S/,[/ii](a:) < 777— 2 ^2(^)1 Mx) := 4 7^(^(1 - z)) 2 + tt^tM 1 " s )) f sh (^(l + 



Consequenty we get 



Ho(x) < 



2! 
D 2 



2-1! 



Step 3. We give an estimate for H m (x) by induction. 
Suppose that for x G [—1, 1], 



# m -i(x) < 



m— 1 



m— 1 



— —-^^(x), /^(z) := sh(z h (l + x)) ^ , 

v ft/ 

(14.251) and this ineqi 
easily see by ffl~26j) and ffl~2Tj) that 



It follows from fl4~25l) and this inequality that H m (x) < b m - 1 %[h m -i]{x)/{2hz 2 h ) m - 1 . We 



5h[h m -i](x) < 



sh(zh(l — x)) d 



2hz h sh(2z h ) 



E a m-l,k 
Om-l-fc^j 



k=l 
m—1 



Sh(z fc (l + x)) v-^ dm-l,* f4( 1_X ) 



Om— 1 



2/iz h sh(2z h ) ^2— i-fe/c! 
Using Lemma I4T5] (14.341) and (j4.35p . we have 
sh(z h (l - x))& k ~ sh(z h (l + x))M k 



\k+l 



k+1 



-ch(2z h ) - M k 



2z 



— <J sh(z h (l + x))(Q 1>k sh(2z h ) + Q 2 , k ch(2z h )) + (Q 2jk ~ Qo,k) sh(z h (l - x)) 



k {z h (l — x)} k ~ l 
By Qi )k + Q2,k = ^2 — «>7T- — ~< — and Q°' k - ^ 2 ' fc ' we get 



1=0 



2\k-l)\ 



bk\ 



k-l 



4<^*W4 W i + ,))£M^. 
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Hence we obtain 



bsh(z h (l+x)) ^ d m . hk j k ^ (z h (l-x)) k - 1 2(z h (l - x)) k+ 1 
b h [n m . l \{x)_ l^ 2m - k yZ^ 2\k-l)\ + (Jfe + 1)! 

Therefore, setting 
h m (x) := sh(^(l + x))^^4y(^(l-^)) fc 

we conclude that 

H m (x) < — — : — h m (x) for x G — 1, 11. 

v ; ~ (2hz h ) m y ' L ' J 

5tep 5. We determine c? TO) fe's. 

It follows from the definition of /i m that for m = 1, 2, . . . , [T//i] and = 2, 3, . . . , m — 2, 

^m,m = dm~l,m-l — ■ • ■ — = 1, G^.l = 1, 

m— 1 k+1 m—1 

(4.37) d m ^ m —\ ^ dl,h d m ,m—k ^ ^j dm—l^m—l) d m \ ^ ^ d m —\\. 

1=1 1=1 1=1 

Here we see that ^ 

cL, m -i =m-l, d m , m -2 = — (m - 2)(m + 1). 

We use these results to obtain d 2 ,2 = 1, ^3,1 = d 3t2 = 2 and d^ 3 = 1. 
By induction we assume that for each m > 3, 

^ fc-i 

d m -i, m -i-fc = — (m - 1 - k) Y[( m - 1 + P) for = 2, 3 . . . , m - 2. 

p=i 

Then we calculate by using (I4.37P and this formula that for k = 3, 4, . . . , m — 2, 

f ' TX\j. 1 r 0\ ^ 

d m ,m-k = l + (m-2) + ^ V — - F! \{{m + q) 

p=2 \P >' q=0 

= (m "T +1) + t j! S^n("'+^ 

p=3 ^ q=0 

{m — 3) -A- (m — 1 — p) P 



3 , n( m +p)+E (p _iv n( m +g) 

p=l p=4 VF g=0 

{m — k)%-^, , (m — fc)(m + — 1)! 
}J_(m + p) 



k\ ' m\k\ 

P =i 
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Replacing k with m — k, we have 

k (2m-k)\ p , 

d m ,k = - : r, Try for k = 2,3, . . . ,m - 2. 

2m — k m\[m — k)\ 

This formula clearly holds for k — l,m — l,m. Thus we obtain (I4.14p . □ 

5 Proofs of main results 



First we prove Theorem I2TT1 Let x$ := minjxi, X4} and 5 G (0, T). Set W\ := [— £5/2, £5/2], 
W2 : = [-^5,^5]) ^6 := m in{L 5 , X5/I6} and 

P fe := [0, M 6 /| log/i|] x Wi, Qs,h := (H x [0, T - 5])\(int P ft U ({0} x Wi)). 

We show that there are i^i, > and h$ > such that for all h G (0, fto)> 

(5.1) sup \u(t,x) - u h (t,x)\ < K 1 Vh, 

(t,x)eP h 

(5.2) sup \u(t,x) -u h {t } x)\ < K 2 Vh\\ogh\. 
Combining these estimates, we obtain the result of Theorem 12.11 

Proof of ( 15. lj) . Choose /i 0jl > so small that h < M 6 /\ log/i| for all ft, G (0, h ,i)- For 
£ G [0, h), Theorem 13.21 (2) directly yields that 

\u(t,x) - u h (t,x)\ = \u(t,x) - tp(x)\ < CVh for all t G [0, h), x G H. 

Hence in the following we consider the case t G Jh ■— [h, Mq/\ logft|]. The u and u h are 
given by, respectively, 

u(t, x) = [T(t)u ](x) — / Ey 5 (t — s, x, x 5 )u(s, x 5 )ds + / Ey 5 (t — s,x, —x 5 )u(s, —x 5 )ds, 

Jo Jo 

m (— ) 

u\t,x) = 9Z,hH(x) + E l(ofl\ ^ sh ^^ - '))](*) 

fc=l s v ^ 5 71 / 



+ E^S)S^[3h(, A (, 5+ .))]W, 

for £ > 0, m = [t/h], x G W2 and h > 0. Here the family {T(t)} t > is a contraction and 
analytic semigroup generated by the operator Au := — u xx + flu in W 2 and D(A) = {«G 
C 2 (W 2 ) I u{±x l ) = 0} (cf. [21, Corollary 3.1.21]). We simply denote E x \ S X5>h by Q h , E, 
respectively if no confusion arises. 

Step 1. We estimate ||[T(£)« ] — S/T[ M o] ||c(w 2 )- 

We use the contraction property of T(t) to have 

||[T(t)tt ] - [«o]||c(w a ) < II [ T (t ~ mh)uo) - «o||c(wy + || [T(h) [u ] - 9h[uo]\\c(w a ) 

+\\[T((m - l)h)[9 h [u )) - gr'lSfcMllccwso 
= : h,i + ^7,2 + 17,3. 
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Since Uq(= (p) is Lipschitz on Q and [T(t — mh)uo\ satisfies u t — u xx + (3u = 0, it follows 
from the theory for parabolic equation that / 71 < Cy/hfoi all t G [0, T), m = [t/h] and 
h > 0. In addition, direct calculations yield that for all ft > 0, 

^7,2 < ||[T(ft> ] -«o||t7(W 3 ) + IM^) - Sft[«o]||o(Wa) ^ ^^Tl. 

As for J 7j 3, we notice that S/i[^o] G -D(^) an d II S/i ['"o] II c 2 (w 2 ) — C/V~h. Since it follows 
from P Theorem 1.3] that ||T((m - l)ft) - < Cft, we get 

i 7 ,s < \\n(m - i)h) - sr^lWShHWcm) < cVh. 

for all m — 1, 2, . . . , [T/ft] and ft, > 0. Thus we obtain 

sup \\[T(t)u ] - % l [u }\\c { w 2 ) < CVh for all ft > 0. 

tE[h,T),m=[t/h] 

Step 2. We estimate I^± : = / -Ej/^ — s,a;, ±xs)u(s, ia^cfts. 
We calculate that 

- s, x, ±x 5 ) - Eo tV (t - s, x, ±x 5 ) | < Ce- x ' /w(t - s) 

for all t, s G (0, T) (i 7^ s) and x G Wi.. It is seen by this estimate that for t G Jh and 
x G Wi, 

|/ 8 ,+| < ||«llc([o,T)xn) / (\E , y (t-s,x,x 5 )\+e- x * /16it - s) )ds 

Jo 

ft 

< M 7 ,i * 5 '^ e-^-^'^ds + M 7i2 h\ \ogh\~ 1 =: J 8>1 + M 7 , 2 ft| logftt 1 . 
./o — S J 

Here and in the sequel the constants M?/s (i > 1) depend on x 5 , but not on ft > 0. 

1 2 2 

Setting r = (x 5 — x)/2y/t — s and using / e~ r cfr < e~ a /a for a > 0, we have 

a 



4,i < M 7t3 / e~ r2 dr < I£^ e -^l/i6t < M? ^ for all te.J h} xeW l and ft > 0. 



Therefore we get J 8j+ < Mf^h for all i G J^, We can get Jg,- < M 7t5 h for all i G J^, 
s G 14 7 ! and small ft > by the same way as above. Hence we have 

\h,±\ < M 7 h for all t G Jh, x G W\ and small ft > 0. 

Step 3. We estimate J 9j± := ^ -^_iL g^-fe[ sh ( Zfe ( Xl T .))]( x ). 

It directly follows from the proof of Theorem 14.51 that |/9,±| < M^ft for all m = 
1,2, . . ., [M 6 /ft| logft|], x G W\ and ft > 0. 
Therefore we have (15 .ip for ft G (0, fto,i)- □ 
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Next we prove (15.21) . The point is to estimate the difference (u m — u m _i)/h — u t . To 
do so, we use the method similar to [UJ, the precise comparison argument of viscosity 
solutions. 

Before proving (15.21) . we recall the definition and some elementary properties of the 
parabolic 2-jets. Let W C M be an open interval. For u : (0, T) x W — > K., we define 
CP 2,± -u(t,x) and CP ' u(t,x) as follows: 



CP 2 ' + u(t,x) 



CP 2 '\(t,x) 




w(t + s, x + /i) < u(t, x) + as + ph + 2~^^ 2 

+o(\s\ + \h\ 2 ) as (s, h) -> (0, 0) 



3{(t n ,x„, ,a n ,p n ,X n )}+r i C (0,T) x W x R 3 such that 

(t n: X n , u(t n , Xn) , Cl n , p n , Xnj > (t, X, u(t, X) , O, p, X^j 

asn-+ +oo and (a n ,p n ,X„) G CP 2 ' + M(t n , x n ) 
CP 2 '-w(t,x) := {(a,p,X) E R 3 \ {-a,-p, -X) E CP 2 ' + (-w(t, x))}, 
9 2, ~u(t,x) :={(a,p,X) el 3 | {-a,-p, -X) E CP 2 ' + (— u(t, x))}. 



We use the following lemma to obtain (15.21) . 

Lemma 5.1 Let u, u t E C((0,T) x W). For any (t,x) E (0,T) x W, if (a,p,X) E 

2 + 2 — 

CP' u(t, x) (or CP ' u(t,x)), then a = u t (t,x). 

Proof. Since u is differentiate with respect to t, we can easily show that for any (t, x) E 
(0,T) x W, if (a,p, X) E CP 2 ' +, u(t, x), then a = u t (t,x). The assertion follows from the 

2 + 

continuity of u t and the definition of CP ' u(t,x). 

2 — 

The case (a, p, X) E CP ' u(t,x) is proved similarly. □ 

Proof of ( 15.21) . First, we show that for any 5 > 0, there exist K 21 > and h 0t2 > such 
that 

(5.3) sup (u(t,x) -u h (t,x)) < K 2 ,iVh\\ogh\ for all h E (0, h Q , 2 ). 



^tep i. We define ;u ft (t, a?) by 
u fc (t,x) := 



Uq(x) for t G [0, h], x E tt, 

u m (x) for t E (mh, (m + l)h], m = 1, 2, . . . , [T//i] and x E Q. 



For any 5 E (0, T), put Tg/ 2 := T — 5/2 and define 



$(t,x,s,y) :=u(t,x)-u (s,y) - —=(t - s) - —=(x - y) 

2\Jh 2vh T 5 /2 - t T 5/2 - s 

Then $ is upper semicontinuous on Qs/2,h x Qs/2,h and $ — ► — oo as t, s S T$/ 2 . Let 
(t,x,s,y) be a maximum point of $ on Qs/2,h x Qs/2,h- We may consider x, s, > 
because if otherwise, we easily get (I5.3P with K 2) i = 4/5. 
Step 2. We show that there is /io,3 > such that 

(5.4) u(t,x)-u h (s,y) < Ch 1/4 \ log h\ 1/2 for all h E (0, h , 3 ). 
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For this purpose, we first study the behavior of (t,x,s,y). It directly follows from 
®(t,x,s,y) > that 

(5.5) -i=(t - s) 2 + — + — ^ < u(t,x) - u h (s,y) < C. 
2Vh T 5 /2 - t T 8 /2 - s 

Hence we get 

(5.6) \t- s\ < Ch l/A . 



Besides, since it is easily seen from $(£, y, s, y) < $(t, x, s, y) and Theorem 13.21 (2) that 

(x - y) 2 < u(t, x) - u(t, y) < C\x - y\ , 



2\fh 
we have 

(5.7) \x-y\<CVh. 

To obtain (15.41) . we divide our consideration into several cases. Set d p Qs/2,h '■= 

dQ6/2,h\({T S /2}_xty-__ 

Case 1. (t,x) or (s,y) G d p Q s / 2 ,h- 

We may assume (t,x) G d p Qs/2,h since the other case can be treated by the same way. 

Subcase 1-1. x G dQ or (t = and |x| > £i/2). 

Then u(t, x) = <p(x) and we get by Theorem 14.11 and (15.71) 



u 



(t, x) — u h (s, y) < <f(x) — (p{y) < C\x — y\ < Cy/h. 



Subcase 1-2. t G [0,Xq/32| \ogh\\ and \x\ < x x /2. 

It follows from Theorem 03] (2), (^Tl) . fl5T6|) and (JEZD that for small h > 0, 

it (£ , x) - u h (s, y) < u (t, x) - u h (t , x) + (t , x) - u h (s, y) 

< CVh + C(h 1/4 ^/\hgh\ + h l/2 ) < Ch 1/A ^/\h^h\. 

Case 2. (t,x), (s,y) G int Q S /2,h- 

Using the maximum principle for semicontinuous functions, we can find a, b, X, Y G 
satisfying 



[a 



,(x-y)/Vh,X)eV 2 ' + u(lx), (b,(x-y)/Vh,Y)eV 2 '~u h (s,y), 



h~\0-Yj-h\-l 1 ; ' (T 5/2 - t) 2 (T 5/2 - s 



Set m = [s/h). Then u h (s,y) = u m (y). In addition, ((a; — y)/Vh,Y) G J 2 ' u m (y) 

2 ± 

and X < K (See [8] or [17] for the defintions of J ' and the maximum principle for 
semicontinuous functions). 

We estimate the difference (u m (y) — u m ^i(y))/h. It follows from $ t (t,x,s,y) = 0, 
( 15. 8p . Theorem 13.31 and Lemma [5.11 with W = Q that 

/- — \ 1 i- \ \h 1 /- _n V~h 

5.9 a = u t {t, x) = —=(t-s) + — =-, b = -={t - s - — 

Vh [Ts/2 -ty Vh {T5/2 - s) 2 
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Noting b > by Theorem 14. 14 we get t > s. We see from ( 13.81) . the first formula of (15.91) 
and this fact that |i — s| < Cwh\ \ogh\. Substituting this into the second formula of 
(E3D, we have l/(T s/2 -s) < Ch-V 4 y/\logh\. Thus ${t,x, s - h,y) < ®(t,x,s,y), (1531) 
and this estimate yield that 

(5.10) Umiv) -^ m - iiv) < ^{(t-is-h^-it-m 

h -l/2 h -l/2 



T5/2 -(s-h) T 5 /2 - s 
1 h 3/2 

By the way, since u is a viscosity subsolution of (I2.2p and 

Urn IS Qb viscosity supersolution 

of (12.41) . we have the following inequalities. 

(5.11) min{a — X + j3u(t,x),u(t,x) — tp(x)} < 0, 

(5.12) min ^Mz^M - y + {3u m (y) , u m (y) - <p(y)\ > 0. 

If u(t, x) — f(x) < in (15.111) . then we easily have by the above inequalities and (I5.7P 

u(t,x) - u m {y) < <p(x) - <p(y) < CVh. 
Thus, in the sequel we assume u(t,x) — (p(x) > for small h > 0. Then by (15. lip . 

(5.13) a - X + /3u{t, x) < for small h > 0. 
On the other hand, we easily get from (I5.10p and (15.121) that 

b + CVh -Y + (3u m {y) > 0. 
Combining (15.81) . (I5.13P with this inequality, we have 

(5.14) u(t,x) - u h (s,y) = u(t,x) - u m (y) < CVh. 

Taking /i ,3 > small enough, we conclude that (15 .4p holds for all h £ (0, h 0t3 ). 
Step 3. We improve (15.41) and establish (15. 3p . 
Substituting (15.41) into (15.51) . we obtain 

(5-15) / 10 := + 7fA~^ < Ch-^y/\^h\. 

-1-5/2 — t J- 5/2 ~ s 

We observe from $(s, x, s, y) < §(t,x,s,y), Theorem 13.21 (2) and this inequality that 

1 



2Vh 



(t-s) 2 < C\\ogh\\t-s\ 
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Hence we get \t-s\< CVh\ log/i|. By using this estimate, we improve the estimates in 
Case 1 as follows. 

(5.16) u{t,x) -u h {s,y) < CVh\logh\ 3/2 if (t,x) or (s,y) G d p Q s/2jh . 
Therefore, (15 .4p can be improved in the following way: 

u(t,x)-u h (s,y) < CVh\ \ogh\ 3 ' 2 . 

Substituting this into (15. 5p again, we get I\q < C\ log h\ 3 ^ 2 . Repeating the above argu- 
ment, we have \t — s\ < Cyjh\ \ogh\ and improve (15.161) as 

u(t,x) -u h (s,y) < CVh\ logh\ if (t,x) or (s,y) G d p Q s / 2 ,h- 

Consequently, we have from (15.141) and this estimate 

x, s, y) < Cy/h\ log h\ for all and h G (0, h 0t3 ). 

Choosing a large K 2> i > C + 4/5, we obtain (15. 3p . 

Next, we prove that for any 5 > 0, there are K 2y2 > and ft,o,4 > such that 

(5.17) sup (u h (t,x) -u(t,x)) < K 2 , 2 Vh\\ogh\ for all h G (0, fr , 4 )- 

(t,a:)eQ«,/, 

^ep ^. Let -u h be defined by ( 12. 5p . For any 5 G (0, T), define 



$(t, x, s, y) := u (t, x) - u(s } y) - — -={t - s) - —j=\x - y 



2Vh ' 2Vh ' T m -t T s/2 -s' 

Let (t,x,s,y) G Qs/2,h x Qs/2,h be a maximum point of $. We may consider &(t,x,s,y) > 
0. Note that Q and (JEZj) hold. 

If (t, a?) or (s, |7) G d p Qs/ 2: h, then we see by similar arguments to those in Case 1 of 
Step 2 that 

u h (t,x) - u(s,y) < Ch 1/4 \ \ogh\ for small h > 0. 

Thus we may assume (t, af), (s, y) G int Qs/ 2t h The maximum principle for semicontinuous 
functions yields a, b, X, Y G M satisfying 

(a,(x-y)/Vh,X) eV 2 ' + u h (t,x), (b, {x -y)/Vh,Y) eV 2 '~u(s,y), 



- V o -y ; - ^ -l i ; ' (t 5/2 (t 5/2 - <o 2 ' 



Let m = [t/7i]. Then note that u h (t,x) = u m (x) and hence ((a; — y)j\fh,X) G J 2 ' u m {x) 
and X < y. 

We estimate (u m (x) — u m -i(x))/h. We see by $ s (t, x, s, y) = and Lemma I5TT1 that 

1 _ a//i 1 _ \fh 

5.18 a = -=(1 -s) + — =-^, b = ut{s, y) = -={t - s) - — — . 
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Note t > s by u t (s, y) > 0. Dividing $(t, x, s, y) > $(t — h, x, s, y) by h, we observe from 
this fact that 

Um{ * ] -^ m ~ 1 ^ > ^L- 2 {(t-sy-((t-h)-sf} 

h -X/2 h -l/2 
+ 



> a 



T5/2 - (t-h) T S /2 - t 

h 3/2 



{Ts,2-tf 



The (IOI . (1531) and the fact t > s yield that l/{T s/2 - t) < Ch'^y/l log h\ for small 
h > 0. Using this estimate, we have 

M m (s)-M m -i(z) > fl _ Cv ^ f or sma ll /,> 0. 

Since the remainder is totally similar to Step 1, we have 

u h (t,x) -u(s,y) < CVh\ \ogh\. 

Thus taking ho,4 > small, we obtain (15.171) . 

Taking K 2 : = max{i^ 2 ,i, -^2,2} + 4/5 and ho,2 '■= niin{/i ,3, ^0,4}, we have the result. □ 

We establish the result of Theorem 12. II by choosing h := min{/i ,i, ^0,2}- 
The proof of Theorem 12.21 is similar to [H] , based on the limit operation of viscosity 
solutions due to [3J H] • 

Proof of Theorem [Ml Recall that x h (t) is given by (12.51) , Define 

x(t) := limsupx ?t (s), x(t) := liminf x h (s). 

Notice x(0) = x(0) = x by Theorem 14.71 

We show x* < x in [0,T). Fix t £ (0, T) and x > x(t). Then there exist sequences 
{h n }t=i an d { m n\n=i sucn that asn-> +00, 

h n ► 0, TYl n h n ^ t, X (^n) *^m n ^ X_(t). 

Since x > x mn for large n e N, we get u hn (s n ,x) = u mn (x) = <f(x). Letting n — > +00, we 
have u(t,x) = <p(x) by Theorem 12.11 and thus x*(t) < x{t). 

To prove x < x* in [0, T), we suppose x(to) > x*(t ) + 6eo for some t 6 (0, T), £ > 
and get a contradiction. By the continuity of x* (cf. (12. 3p ). there exists 5 > such that 

(5.19) x(t ) > x*(t) + 5e for all te (t -55, t + 55). 
Choose {/i n }n=L an d { m n}t=i satisfying 

(5.20) h n — > 0, m n h n — > t , x hn (m n h n ) — > x(t ) asn^ +00. 
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Take n G N such that 

\^nh n — t 1 < 5, \x hn (m n h n ) — x(t )\ < e for all n > no- 
Using ()2.3|) . Theorem 14.71 and these facts, we observe that 

x hn (t) > x hn (m n h n ) > x*(t) + Ae for all t > m n h n and n > n . 
This implies that for all n > no, 

(5.21) u hn >ip in Q := (t + 25,t + 45) x (x*{t ) +e ,x*{t ) + 3e ). 

On the other hand, we notice u = (f in Q. 

Fix y G (x*(t ) + 2e ,x*(t ) + 4e ). We derive 

(5.22) -fe(!to)+MW))<0' 

Let = (f>(t, y) be a smooth function such that u — <fi takes its strict maximum at (to + 
2S/5,y ) in Q and (<f>(t + 36, y ), <fr xx (t + 35, y )) = ((p(y ), ip xx (y )). Let (t n ,y n ) be a 
maximum point of u hn — (f> in Q. Then it can be observed from Theorem 12.11 that 

(5.23) (t n ,y n ) — ► (t + 3<y,y ), u hn {t n ,y n ) — ► u(t + 35, y ) = ip(y ) as n -> oo. 

Put m n = [i n //i n ]. Then from (15.211). « in (|/ n ) = u hn (t n ,y n ) > ip(y n ) for n > n . Using 
the fact that is a viscosity subsolution of (12. 4p with m = rh n and Theorem 14.11 we 
have the following inequality. 

-fxxiVn) + f3um n {y n ) < 0. 

Letting n — * +oo, we get (I5.22p . 

However, (15.221) contradicts to (13.101) because of yo > x*(t ) > xq. Thus we have 
x < x* in [0, T) and conclude that x = x = x* in [0,T). Applying jSJ Section 6], we 
complete the proof. □ 

6 Appendix 

6.1 Formal asymptotic expansion for (12.41) 

This subsection is devoted to the formal asymptotic expansion of the solution of (12.41) 
with m — 1 near the free boundary as h \ 0. 

Let x* > be the free boundary of (12 .4p with m — 1. From the facts U\ = tp on [x*, 1] 
and U\ G C 1 (fi), it is sufficient to treat the following problem instead of (J2T4]) : 

(6.1) U - - - u xx + (3u = for x < x*, u(x*) = y?(x*), u x (x*) = ^(x*), 

We rewrite (I6.ip . Set w(x) := (w(x) — ip(x))/e ax . Then w satisfies 

(6.2) —hw xx — 2ahw x + (1 + r/i)w = h(—qe x + r) for x < x*, w(x*) = w x (x*) = 0. 
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The solution w of this problem is given by 

(6.3) w(x) = k ie x+{x - x * ] + k 2 e x - {x ~ xt) + h (ZS^L + —L—] for some fc ls k 2 E R, 

\1 + qh l + rhj 

where \± := —a±z h . Since we see from ( 16. 2 p that w(x) = 0(h) as h — > 0, we have k 2 = 0. 
Moreover, we observe by the conditions for u at x* and Taylor expansion to log(l + s) 
around s = that as h \ 0, 

, R/1 x * / (z h -a)(l + qh) \ ( Vh 

(6.4) x = x + log < — — ■ > = x + log 1 + 



(z h - a - + rh) J to ^ v/l + Z^ + av 7 /^ 

*o + ^ =-if ^ _] +0(^/2) 

VI + « v 7 ^ 2 I Vl + ^ + 



+ Vh- (a + l/2)h + 0{h 3 ' 2 ), 



h = ^— + 0(h 2 ). 
l + qh v ' 

where xq is given by (12.31) . Therefore using (16.31) with these results, we have u(x) —<f(x) = 
0(h 3 / 2 ) near x < x* . 

To obtain the asymptotic expansion of w in terms of h, from the above estimate, we 
may assume that w can be expanded as follows: 

(6.5) w{x) = h 3 l 2 w 3 (p) + h 2 w 4 {p) + 0(/i 5/2 ). 

Here p := (x — x *)/2Vh. We impose the following from the conditions for w at x* in (16.21) : 

(6.6) Wi{0) = io-(0) =0 for i = 3, 4. 

Substituting (16.51) into (16.21) . we have by (3 = a 2 + r and Taylor expansion to exp(x + 
Vh - (a + l/2)h + 0(h 3 / 2 )) as /i \ 

(6.7) h 3 ' 2 f-^- + w 3 \ + h 2 f-^ - aw' 3 + u> 4 ) 

= -r{h 3 / 2 {2p + 1) + h 2 (2p 2 + 2p-a)} + 0(/i 5/2 ). 

Here we have neglected kie x+y since this is smaller than the last term of (16. 3p for x < x* 
and h \ 0. 

We determine w^s (i = 3, 4) from the above expansion. Comparing both sides of (16.71) . 
we can derive the following: 

w" w" 
w 3 -^- + r(2p + 1) = 0, w 4 -^--aw' 3 + r (2p 2 + 2p-a)=0. 

Solving these equations under (16.61) we obtain 

w 3 (j>) = r(e 2f} - 1 - 2p), w 4 (p) = r{e 2 " - (1 + 2p + 2p 2 ) + a(e 2p - 2pe 2p - 1)}. 
Therefore we conclude that as h — > 0, 

u(x) = ip(x) + e ax \h 3 ' 2 w 3 (^=P) + h 2 w A (^f~~) + °^ 5/2 ) } near x < x *' 
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6.2 Proof of ( jurrp 

First, we may assume m > 3 and consider the case k — 1, 2, . . . , m — 1. Because the case 
m = 1, 2 or A; = 0, m is easily proved. We see from Stirling's formula that for all p € N, 

+ I2p~- v/27rpP+i/2 e -P - + 12^ + ^ 2 ' 



where C is independent of p. Using this inequality with p = m, 2m — k,m — k, we observe 
that for k = 1, 2, . . . , m — 1, 

(2m -fc)! < 1 + C/m { m \ " (m+1/2) f m-k \ ^ m ~ k+l ^ 



2 2m - k m\(m - k)\ ~ ^ (2m - k)n \ m ~ k / 2 J \m-k/2 



Set 5 = 8(m, k) := kj \j2m — k. Then we have 



\ m+l/2 / , \ m-fc+1/2 / c- \ m+1/2 / r \ m-fc+1/2 

m \ I m — k \ I _ o x x A x 



1+ . 1 



m-k/2) \m-k/2) \ y/2m - k J \ y/2m - k. 

Denote the RHS of this formula by In. Since it follows from Taylor's expansion that 



log(l + x) = x - ^x 2 + ^x 3 - 4( - 1 ^^ 4 , |0| < |z| < 1, 

we obtain for k = 1, 2, . . . , m — 1, 

. 1W 5 5 2 1 5 3 1 rt : 

login > ( m+ - ] ( - — — + 



2 J \ y/2m - k 2(2m-fc) 3 (2m - /c) 3 / 2 4 (2m - /c) 2 

, 1\ ( 5 5 2 1 5 3 1 5 4 

+ m — k H — 



2y V V2m- k 2(2m-k) 3 (2m - A;) 3 / 2 4 (2m - A;) 2 

1.2 ^ <^ <^ / 1 

:5 - ^ TT + 7T ~ T7^ 7T 1 + 



2 2(2m - k) 3(2m - k) 4(2m - k) \ 2m - k J ' 

Since it is easily seen by m > 3 that 

11/ 1 \ 1 1 

1 + > > 0. 

3 4 V 2m - k ) ~ 12 4m ~ 



Thus we get 

1 ^ 5 2 1 



log In > - — > -5 2 . 



2 2(2m - Jfc) 4 
Therefore we obtain (14.171) . □ 
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